Super-rough phase of the random-phase sine-Gordon modeh Two-loop results 
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We consider the two-dimensional random-phase sine-Gordon and study the vicinity of its glass 
transition temperature Tc, in an expansion in small r — {Tc — T)/Tc, where T denotes the 
temperature. We derive renormalization group equations in cubic order in the anharmonicity, 
and show that they contain two universal invariants. Using them we obtain that the correla- 
tion function in the super-rough phase for temperature T < Tc behaves at large distances as 
{[e{x)-e{0)]^) = Aln^{\x\/a) + 0[ln{\x\/a)], where the amplitude A is a universal function of 
temperature A = 2t^ — 2r^ -I- ©(r^). This result differs at two-loop order, i.e., C'(r"^), from the 
prediction based on results from the "nearly conformal" field theory of a related fermion model. We 
also obtain the correction-to-scaling exponent. 

PACS numbers: 64.70. Q-,64.60ae 



I. INTRODUCTION 

Although in two-dimensional (2D) systems with con- 
tinuous symmetry and short-range interactions thermal 
fluctuations prevent the existence of long-range order—, 
they do not prevent phase transitions. The 2D XY 
model, much studied in that context, describes a large 
class of physical systems with continuous symmetry, 
which includes superfluid and superconducting films, 
magnetic systems and one-dimensional quantum liquids. 
It exhibits a topological phase transition between the 
low-temperature phase with quasi-long-range order and a 
disordered phase at high temperatures. This Berezinskii- 
Kosterlitz-Thouless transition is driven by unbinding 
of vortices due to an increasing amount of thermal 
fluctuations^!^. From the technical side, the XY model 
is conveniently studied within the equivalent dual 2D 
sine-Gordon (SG) model, which is amenable to powerful 
field-theoretical treatments^. It exhibits a high temper- 
ature quasi-long-range ordered phase and a low temper- 
ature massive phase. When additional terms are added 
to the SG model new universality classes can emerge^. 
For a pure system, the simplest example is an addi- 
tional field gradient in one direction^, which describes 
the commensurate-incommensurate transition in 2D and 
realizes, for example, when an atomic layer of noble gases 
is deposited on the periodic substrate of graphite^. Both 
models, with and without the tilt, are exactly solvable^ 
and are by now well understood. 

The random versions of the 2D SG model allow for 
more scenarios and much less exact results are known. 
Via bosonization they are related to fermions with disor- 
der and have also been much studied in that context^^. A 
well known example is the 2D SG model with a quenched 
random phase that depends on only one coordinate in the 
cosine term. Such a model describes a classical 2D model 
with correlated disorder, and also a ID quantum system 
with point disorder (the second direction being imagi- 
nary time). In the latter case it is related to ID disor- 
dered Luttinger liquids and belongs to the Berezinskii- 
Kosterlitz-Thouless classic. Seen as a 2D classical model 



it exhibits quasi-long-range order (i.e., an infinite corre- 
lation length) in its high-temperature phase, which is de- 
scribed by a line of Gaussian fixed points of the renormal- 
ization group (RG) where the cosine term is irrelevant. 
Its low-temperature phase is glassy and described by RG 
fixed points at large disorder strengths, see Fig. [21 This 
scenario first found in one-loop order— is not changed at 
two-loop order—. 

In this paper we study the 2D random-phase sine- 
Gordon (RPSG) modeliiii^ii^ where the quenched ran- 
dom phase depends on both coordinates. This model 
can also be understood as the random field XY model 
provided one excludes vortices by handi^. The RPSG 
model describes 2D periodic disordered elastic systems, 
such as a randomly pinned planar array of vortex linesi^ 
or surfaces of crystals with quenched disorder ^i. It 
also exhibits a phase transition at a critical tempera- 
ture Tc below which the random cosine term becomes 
relevant. This transition was first studied in an expan- 
sion in r = {Tc — T)/T using a one- loop RG approach 
in Ref.^"'. The physics of the RPSG model is however 
quite different from the previously mentioned disordered 
model. While the high-temperature phase is described 
by a line of Gaussian fixed points, similarly to the SG 
model, the low-temperature phase is glassy and described 
by a line of non-Gaussian fixed points where the renor- 
malized disorder gradually increases from zero when de- 
creasing the temperature below T^ see Fig. [T] The glass 
phase for T < Tc is super- rough, i.e., the variance of the 
fiuctuations of the displacement field grows as logarithm- 
squared of the distance in contrast to the standard rough 
logarithmic form at high temperaturesii. Within the 
RG approach this is due to an unbounded growth of 
the off-diagonal disorder (in replica space) at low tem- 
peratures that determines the correlation function. Fur- 
ther RG studies predicted that the amplitude A of the 
logarithm-squared correlations is a universal function of 
temperatursii^, with A — 2t^ -|- 0{t^) to one loop order 
accuracjJii^. 

The existence of the super-rough phase has been 
confirmed in several numerical studies at zero 
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FIG. 1. The renormalization group flow diagram of the 
random-phase sine- Gordon model U]). The line of flxed points 
(red curve) in the low-temperature phase T < Tc occurs at 
flnite disorder strengtli and it is continuously connected with 
the line (blue curve) of Gaussian fixed points at T > Tc. 
While the line of nonzero fixed points is a linear function of 
temperature to the lowest one loop order, it gets quadratic 
correction beyond that. The correlation function at T < Tc 
has a super-rough logarithm-squared formii. 
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FIG. 2. The renormalization group flow diagrams of the sine- 
Gordon model. At small temperatures, T < Tc the tentative 
fixed point is disconnected from the line of Gaussian fixed 
points that exist for T > Tc when the strength of the anhar- 
monic term g of the SG model flows to zero. In the context of 
one-dimensional fermions with disorder the same flow diagram 
apply provided the temperature is replaced by the strength 
of interaction and that g denotes the disorder strengtbiS. 



temperatur o^^'^" and for all temperatures < T < Tc, 
see Refs.— These studies consider discrete random 
height models, or discrete-line models, believed to be in 
the same universality class as the RPSG model. These 
are further mapped onto the dimer covering problem 
with random weights^"—. Powerful polynomial algo- 
rithms then allow to generate all possible coverings of 
the lattice by dimers^S.. Using such algorithms A{t) was 
estimated in^i where the quadratic behavior A(t) cx 
at small r was confirmed. Very recently more accurate 
data have been obtained iv^^^ (see below) . 

Some recent studies opened the hope that A could 
be obtained non-perturbatively. Considering a model of 
disordered noninteracting fermions in 2D, Guruswamy, 
LeClair, and Ludwig^* used methods of "nearly confor- 
mal" field theory to predict the exact form of the corre- 
lation functions as well as the scaling equations for their 
fermionic model. Upon bosonization, these results where 
interpreted as corresponding to the RPSG model exactly 
on its fixed points at finite disorder (see Fig.[T]). In Ref.^^ 
this correspondence and the translation to the parame- 
ters of the RPSG model was performed in details, with 
the conclusion that if the exact beta function of^^ is cor- 
rect then one should have A = Ancft = 2t^(1 — r)^ 
exactly in the whole super- rough phase < t < 1. As 
discussed in^^ this however raises some puzzle: numerics 
exclude the amplitude vanishing at T = 0, and the non- 
monotonous behavior of Anc ft with temperature is sur- 
prising. Hence the formula, correct to one loop accuracy, 
can hold exactly at best in a vicinity of T^ i.e., for t < t* 
with some unknown r* . Since the numerical values are 
larger by a factor « 4 than the maximum ANCFTi^/"^) 
the true amplitude should be larger than the predicted 
one. Possible scenarios are discussed in Refi^ such as 
the mapping between free fermion models and the RPSG 
model failing below some temperature, or some new op- 
erators becoming relevant at r* . In addition a functional 
RG study performed in Ref.?^ leads to a non-vanishing 



amplitude at T = as a result of including higher har- 
monics of the disorder that are relevant there. 

In the present paper we revisit the model using per- 
turbative renormalization group methods to the next 
two-loop order and compute the amplitude of correla- 
tions in the super-rough phase including 0{t^) terms. A 
short summary of the present work has been presented 
in Ref.^. Here we give all the details. We perform a 
systematic calculation in terms of the strength of the an- 
harmonic term g, see Eq. ((8]). The problem is studied 
within a bosonic formulation using field theory methods. 
We use two complementary methods, which are explained 
in a pedagogical way. The first one is based on the cal- 
culation of the effective action and the second on the 
operator product expansion. We study both the theory 
regularized by a small distance cutoff a and, for T < 
directly in the continuum limit, and obtain the precise 
dependence of the results on the cutoff functions. Our 
main findings are the scaling equations (I76|) and (j77|) . the 
correlation function (I85p and the correction-to-scaling ex- 
ponent ([50)1 . The equations beyond lowest order contain 
non-universal coefficients that are connected by relation 
([75)1 . For the amplitude of correlations in the glass phase 
T < we find A = 2t'^- 2t^ + 0{t'^). Hence it confirms 
the conclusion o^ that the translation of the results of 
Ref.^'^ into an exact result to all orders for the RPSG 
model cannot be correct. Since the discrepancy arises 
already at two loop order (i.e., r* = 0) it also casts some 
doubts on the fermion calculation done in Refi^ or on 
its consequences for the bosonic model as given in Ref.— 
and in Ref.^-. Note that our result for the amplitude is 
indeed larger than Ancft (for small r) hence it goes 
in the right direction. Interestingly, it also appears to 
better fit the most recent numerical results up to 

T w 0.5, although it is a perturbative calculation around 
Tc, i.e., around r — 0. 

To be complete let us mention that many other 
studies have addressed the RPSG model, its 
thermodynamic o^^'^" , its stability to RSB and links 
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to fermions^, its dynamics near T,r^—: its equilibrium 
dynamics at all T— and its aging dynamicsiSi3Si2S. It 
would be interesting to push such methods to two loop 
accuracy, as done here for the statics. 

The outline of this paper is as follows. In section II we 
introduce the model. Using the replica method we derive 
the replicated Hamiltonian which is our starting point for 
the systematic field-theoretic renormalization group pro- 
cedure. We also define several correlation functions of 
interest. In section III we calculate the effective action 
of the model order by order in the disorder strength and 
from it we derive the beta functions. Further we examine 
the universality of coefficients in the beta functions. In 
section IV we evaluate the coefficients in the beta func- 
tions using two different methods and find three universal 
coefficients, one of them from one-loop and two of them 
from two-loop. In section V we give the final form of 
the scaling equations, find the correction-to-scaling expo- 
nent and obtain the correlation function (that measures 
the fluctuations due to the disorder) in the super-rough 
low-temperature phase. In section VI we present a first 
principle derivation of the correlation function. In sec- 
tion VII we use the operator product expansion method 
which is found in agreement with the results of the effec- 
tive action method. Section VIII contains conclusions. 
Numerous technical details are relegated to appendices. 



II. MODEL, REPLICATED HAMILTONIAN, 
AND CORRELATION FUNCTIONS 



A. Model 

We consider the 2D random-phase sine-Gordon model. 
In terms of a real displacement field 0{x) G (—00, 00) its 
Hamiltonian reads 



H 



1 



h-VJ~- {ie'^ + h.c.) 



(1) 



averages by . . .. Depending on context x and y will 
be used either to denote 2D coordinates (as in previous 
equations) or as their norms, i.e., x stands either for x 
or Ixl. 



B. Replicated Hamiltonian 

We use the replica method to treat the disorder 
The partition function for the model ([1]) is given by 
Z = / 'D6e~^/'^ . In order to perform the disorder av- 
erage, we use the replica trick, and the free energy of the 
system = — Tin Z is written as 



F = lim (Z" 

n-i-O 



l)/n. 



(5) 



The average with respect to disorder can now be done 
since one can write Z" — J (00=1^^")^^^ ''^^ where 
Oa,a — 1 . . .n are the replicated fields. In the following 
by greek indices a, /3, . . . we denote replicated fields and 
we do not write explicitly the boundaries in the sums. 
The replicated Hamiltonian reads 



where the harmonic part is 



rep 



(6) 



TTTep 

-"0 
T 



E / d'^{^^"^ [(V.0a)2+m2(0„) 

Q/3 



Ait 



(7) 



The mass m is introduced in the model as an infrared 
cutoff. We will perform calculations with finite m and 
study the limit m — at the end. The system size 
is infinite throughout the paper. The anharmonic part 
reads 



rep 



/ d^a;cos( 



a/3 



0- 



(8) 



where n is the elastic constant, a is the short-length-scale 
cutoff, and h{x) and ^(x) are quenched Gaussian random 
fields, the first one real and the other complex. Their 
nonzero correlations are given by 



W{x)ho{y)^ ^h5''5{x-y), 



ax)C*iy) = A^Six-y), 



(2) 
(3) 



where i,j E {1,2} denote the components of h. For fu- 
ture convenience we define the disorder strengths A^^^ in 
terms of the dimensionless parameters a and g as follows: 



27r' 



2 9_ 

2^' 



(4) 



We introduced the symbol ^ which denotes a summa- 
tion where all replica indices are different. While after 
replicating the model ([l} one formally obtains a sum over 
all unconstrained replica indices, for convenience we use 



Correlation functions 



Our aim is to compute the two correlation functions: 



g{x) ^ {e{x)e{Q)) - {e{x)){em, 



g,,{x) = {e{x)){em, 



(9) 
(10) 



where T is the temperature. Note that the disorder h{x) 
must be introduced as it is generated by the symmetry- 
breaking field under coarse graining. We denote disorder 



where Q{x) measures the (disorder averaged) thermal 
fluctuation while ^^0(2;) measures the fluctuations due to 
disorder of the (thermally averaged) displacement field. 
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These disorder averaged correlations can be obtained 
from correlation functions of replicated fields. For in- 
stance 



g{x) = lim - y^Gapix), 



lini(l - 6al3)Qaf}(x)., 
n— >0 



where 



It can also be expressed as 



(11) 
(12) 

(13) 
(14) 

which contains both correlations defined in (|9]). G{x) is 
also called the connected part and Qq [x) the off-diagonal 
part. To this aim we will use the harmonic part Hq'^^ as 
the " free" theory and treat in perturbation theory, 
i.e perform a perturbation theory in g. Here and below 
we denote by ((..)) (exact) averages over the complete 
Hamiltonian H'^'^p and by (..) averages over the free part 
Hq'^^ (it also designates thermal averages in the unrepli- 
cated theory, as no ambiguity can arise). 

We start by computing the correlation function for the 
harmonic part, i.e., for g — Q. It is easily found in Fourier 
space: 



T 



1 



271 {q^ + m?) 



2\2 



0{n), 



(15) 



where in the last step we inverted a replica matrix keep- 
ing only non- vanishing terms in the replica limit n — >■ 0. 
Going to real space one obtains 



Go.p{x) = {e^{x)ep{Q)) - 5^,iG{x) + Go{x). 
The connected part reads 
G{x) 



(16) 



2(l-r)ii:o(mVa;2+a2), (17) 

and everywhere in the paper the parameter r denotes 

T = 1 - T/T,. (18) 

By K[) we denote the modified Bessel function of the 
second kind^. Our model ([ij has a phase transition 
at temperature = ^ttk (see Fig. [T|), i.e., for t — Q. 
In the following we will repeatedly use the expression for 
G{x), which has the following behavior at small distances 
|a;| ^ (cm)^-'- 



G{x) = -(1 - t) In [c^m^ix^ + a^)] , 



(19) 



with the constant c = e^^/2 and is the Euler con- 
stant. In (IT71) we have introduced the ultraviolet regular- 
ization by the parameter a. Such choice of regularization 



is preferable to some other choices in momentum space, 
since our RG procedure is most easily done in coordinate 

41 

space.— 

The off-diagonal part of the correlation function for 
Q = reads 



Go(x) 



(27r) 



K(3{m\x\) 



^-Ki{m\x\) 



0(n), 
(20) 



which at small distances |a;| ^ (cm) ^ becomes 

Go{x) = -2cr(l - rf {1 + In [c^m^{x^ + a^)] } + 0{n). 

(21) 

The model studied here possesses an important sym- 
metry, the statistical tilt symmetry (STS), i.e., the non- 
linear part H^'^ is invariant under the change O^ix) — >■ 
da{x) + (j){x) for an arbitrary function (l){x). As discussed 
in many worksi^iii^i^, and recalled below, this implies 
two important properties: 

(i) Go (x) does not appear to any order in perturbation 
theory in g in the calculation of e.g. the effective action 
(see the following section) , and (ii) the disorder averaged 
thermal correlation is uncorrected to all orders: 



g{x) = G{x), 



(22) 



i.e., independent of g. This implies that t or T/Tc can be 
measured from the amplitude of the logarithm in g{x) ~ 
2(1 — r)lna; at large x, hence they are uncorrected by 
disorder 4^ 

Because of property (i) Go{x) only receives addi- 
tive corrections, e.g. corrections to a which, in the 
present model, change its logarithmic behavior (l?T|) into 
a squared- logarithm behavior for Qq^x), as discussed be- 
low. 

The perturbation theory thus depends only on the 
function G{x). While the precise form of the correlation 
function for the model ([7]) can be explicitly calculated 
and takes the form (|17p . we will explicitly check below 
that the precise form of G{x) is not important as long as 
it satisfies the two conditions: (a) the limiting behavior 
of the correlation function (|17p is logarithmic as given 
in ([T9|) and (b) the propagator tends (exponentially) to 
zero at large distances. The crossover length is given by 
the infrared cutoff, which is the inverse mass in our case. 
The freedom of the propagator that satisfies conditions 
(a) and (b) is manifested through the renormalization 
group equations which will contain several non-universal 
constants that, when appropriately combined, produce 
some universal numbers. These universal numbers deter- 
mine in turn the amplitude of the correlation function 
in the super-rough phase, as well as the correction-to- 
scaling exponent. 



III. EFFECTIVE ACTION OF THE MODEL 

In this section we calculate the effective action func- 
tional r for the model ([T]). It will directly lead to the 
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scaling equations of the model ([T]) and critical properties 
of the system. It extends to the next order in perturba- 
tion theory the calculation of F in^^. In the framework of 
diagrammatics, F can be expressed as a sum over all one- 
particle irreducible graphs^'^. Here we will not calculate 
F using a diagrammatic approach but via an equivalent 
algebraic method. The definition and the derivation of 
the final form of the effective action (|A17I) . to the re- 
quired order in perturbation theory, is presented in ap- 
pendix El The difference with the standard Wilsonian 
procedure of Ref.— is that one integrates out fields that 
live in the whole momentum space and not only in the 
high-momentum degrees of freedom. These fields that 
are integrated out are denoted by x in (jAl7p . Translated 
to the replicated Hamiltonian, our aim is to evaluate the 
following expression 



TTTep 

-"o 
T 



Fi+F2+F3 + 0(/), 



(23) 



where F^ is the corresponding term from (|A17|) propor- 
tional to g^. 



Derivation of F 



To lowest order in g we have 
T, = {Hl'^{e + x)/T)'' 



27ra2 



(24) 



Q/3 



In the previous and all forthcoming terms of similar form, 
evaluation of averages of the type (• • • is quite simply 
done by making use of Wick's theorem when one gets 
contractions of x fields with respect to the quadratic 
Hamiltonian ([7]). Due to "charge neutrality" of ([8]) only 
the diagonal part p7)) of the correlation function Gq^ (x) 
survives in expressions of the type (^[Hl'^^{9 + x)/T]^)'^, 
where p is a positive integer. 

In order to obtain F2 we use the transformation of sum 
rule (jBip when evaluating the corresponding term from 
(|A17p . The final result reads 



F2 = " - 



27ra2 



-2G(0) 



s=±l (. aP 



J2 \J2 d^xd^yA{x-y,2s)cos[9c,{x)-se^{y)-0p{x)+sep{y)] 



where 



J2 I d''xd'y2Aix~y,s)cos[ea{x)~sea,iy)-ep{x)+se^{y)] }, 



A(x,p) =eP^(^) -l~pG{x). 



(25) 



(26) 



The second-order term F2 consists of two- and three-replica contributions. It turns out that only terms with s = +1 
give contributions to the renormalization. The term multiplied by A{x — y, 1) is responsible for renormalization of 
the coupling constant g, while the other term multiplied by A{x — y, 2) renormalizes the off-diagonal part of ([7]). 

One should notice that the second term of Eq. (|A17p proportional to g^ that contains the integral basically makes 
the total contribution (proportional to g"^) of Eq. (jA17p to be one-particle irreducible and produces terms —aG{x) in 
Eq. (Uni)- 

In order to obtain the next-order contribution, similarly as we did for F2, we use the sum transformation (jB2p when 
we evaluate the term (jA17p proportional to 5^. The final result reads F3 = F3 + F3 with 



F' 



-3G(0) 



d?xd?yd? z 



^ 6 V 27ra2 , 

iB{x -y,y-z,z-x, -2, 2, 2) ^ ' cos[0„(a;) -I- ^^(y) - 0„(z) - 9p{x) ~ 9^{y) + 9piz)] 

afi 

+ l2B{x - y, y - z, z - 2, 1, -1) ^ ' cos[9a,{x) - 9c,{y) + 9c,{z) - 9p{x) + 9p{y) - 9^{z)] 
+ 2B{x -y,y- z,z-x, 1,1, cos[6l„(x) - 9a{y) - 9p[x) + 9fi{z) + 9^{y) - 9.y{z)] 
+ 6B{x -y,y- z,z-x, 1,1,0) J2' cos[6i„(x) 9a{y) + 9p{y) - 9p{z) + 9^{z) - 9s{x)] I, 



(27) 
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while the other part Fg is not important for critical properties of the model and is given in (jCl[) . In the previous 
equation the common term B is defined as 

B{x, y, z, a, b, c) =e'^G(.)+6G(a)+cG(^) _ ^aOix) _ ^bG{y) _ ^cG(z) + 2 - Bi (x, y, z, a, b, c), (28) 

where the part of B that makes it one-particle irreducible reads 

Bi {x, y, z, a, b, c) =e'^G(x) j^^^^^^ ^ ^^^^^j _^ ^bG(a) [^(^(^) ^ ^g-^^^j _^ ^cG{z) j^^.^^^ _^ ^^^^^j 

- abG{x)G{y) - acG{x)G{z) - bcG{y)G{z) - 2aG{x) - 2bG{y) - 2cG{z). (29) 



It arises from the last two terms in Eq. (|A17p . 

The summands of Fa can be distinguished by the sum 
a+b+c of the corresponding _B-functions. For purposes of 
renormalization of the model ([T|) the only relevant terms 
are these for which a + b + c equals either two or three. 
The former contribute to the renormalization of the cou- 
pling constant g, see Eq. (|8]), while the latter renormalize 
the off-diagonal part of ([7]). All other summands pro- 
duce nondivergent contributions to the effective action 
and hence can be neglected. For completeness of presen- 
tation they are given in Eq. (jCl[) . 



two-replica part of F2 [Eq. ((25|) ] one obtains 



B. Expansion of T 

Having obtained the effective action F in the preceding 
part, we are now prepared to study its renormalization. 
The perturbative expansion of F in the bare coupling 
constant g contains divergencies when the cutoff a tends 
to zero. In order to remove such divergencies from the 
theory it turns out that two renormalization constants 
suffice. They relate the initial coupling constants a and 
g and the "renormalized" ones, an and gu respectively. 
Written in terms of renormalized quantities, the effec- 
tive action will be free of divergencies (up to third order, 
which is the order we are working with) . In the following 
we will calculate the divergencies of the effective action 
in a double expansion in two small parameters, g 
and r, meaning close to the critical temperature and for 
weak anharmonic terms of the model ([T]). 

We now write F in the same form as the starting Hamil- 
tonian with new coefficients, plus irrelevant terms. That 
allows us to define gn and an below. 

The first-order term (|24p is already in a proper form. 
The second-order term F2 [Eq. ((25|) ] contains terms with 
two and three replica sums. The former gives contribu- 
tion to the off-diagonal part of Hq'^^, while the latter 
changes H\^^ , as we will see below. It is important to 
stress here, that even if we had started without the term 
^ h in ([1]) [i.e., without the term ^ cr in ([7])] this term 
would have been generated under the RG coarse grain- 
ing procedure, as first noted by Cardy and Ostlund— . 
This term is very important for the behavior of correla- 
tion function at T < Tc and its super- rough In^ form.^^ 
We will come to that point later when we investigate the 
correlation function of the model. After expanding the 



9 



'2 [ 27ra2 
27ra2 



-2G(0) 



E' 

Q/3 



J d'xiv^ea^-v^epf 



d^xcos{9a — Ojs) 



(30) 



where . . . stands for many irrelevant operators and we 
define the dimensionless integrals: 



ai 



a2 



9 ? 
cm 

2tt 



2tt 



dV(y,i), 

d^yy^A{y,2). 



(31) 
(32) 



Finally the cubic term in g after expansion produces 
the following terms 



2 V 27ra2/ 



-3G(0) 



61 - 862 + 12af / 2 



a/3 

d''xcos{9a — dp) 



■ j d?x{^^eo.~^jpf + ...^ 



where 



4 4/" 

bi J d^xd^yB{x + y, x, y, -2, 2, 2), 

4 4/* 
C 771 I 

^^"(2^ J d^2;d2yB(x,2/,a;-^j/,2,l,-l), 
J d'^xd^yx'^B{x,x + y,y, 1,1,1). 



b.^ 



(2^) 



(33) 

(34) 
(35) 
(36) 



are dimensionless integrals. We emphasize here that in 
the above expressions (|30p and (p3|) we already take into 
account the replica limit n 0. Had we kept n we would 
have obtained 2(2 — n)ai instead of 4ai in ([50)) and 61 -f 
A{n-2)b2 + 2{n-2){n-3)al instead of &i - 862 12a^ in 
([33| . We should also mention that the term in (|30|) is 
basically the second-order term J d^xd^yB{x, y, 0, 1, 1, 0) 
that could be written as repetition of (PT|) term. 
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The final expression for tlie effective action reads 

~ ^y.Oo, ■ y.9p ~ ^c^TO^ cos(0„ - 9p)} (37) 

47r ZTT i 



where 



gn ^g~2a{g' + - 8&2 + 12a2)53, (38) 

02-2 I ~3 



(39) 



with g = ge "-^'^"^/(cma)^. In the following it will be 
useful to have the inverse expression of ([55]) which reads 

g = 5fl + 2ai5^ - - 862 - ^aDg^. (40) 

In p7p we have returned to the unrestricted replica sum, 
since Y.'^fi{V Jo.-^ xOpf = Y^o^pi^ xOo.-^ xOpf , while 
the cosine term ([5]) differs from the corresponding cosine 
in ([57]) by ©(n) which goes to zero. 



C. Beta functions of the model 

In this subsection we obtain the general form of the 
beta function of the model ([T]) in terms of the integrals 
a,,6i defined in and ([Ml) -(ESI)- We are gener- 

ally interested in the fiow of the effective action when the 
cutoff is varied, and the beta functions describe the flow 
of the terms (i.e., operators) which become relevant at 
the transition. Here they are thus defined by computing 
derivatives of gR and gr with respect to the cutoff for a 
fixed microscopic model, i.e., keeping g and a fixed, as a 
function of gn itself: 



- md^gR = PgigR), 

- mdmO-R = Paign)- 

We thus obtain from ([5^1) : 
Paign) ^'2TgR - 2(2aiT - mdmai)gR 



(41) 
(42) 



+ _[4r(6i -8&2+4a2) 

- mdmih - 8&2 + 4.al)]g% + 0{g%), (43) 

Pa{gR) = ^2Ta2 - ^mc),„a2^ g\ 

+ [2r(4aia2 - 363) - 2aimdyna2 + ^5^63] g\ 
+ 0{g\) (44) 

When deriving the last two expressions we used that in 
the limit ma — )• one can replace^ = g{cma)~^'^ which 
gives mdmg = —2Tg. We have also used the inverse 
relation (|iO|) . 

Eqs. (j331) and (|44p should be understood as an expan- 
sion of beta functions in powers of gR where the first 



three powers are taken into account. Anticipating that 
the fixed point is gR ^ 0{t) for r > we thus need to 
compute the coefficients of gjf to 0{t) and the ones of 
gfj to order 0{t^) to study this equation in the vicinity 
of the fixed point consistently to the desired order in r. 



D. Universality and the beta function 

From the above considerations we can thus surmise, 
and will check below by explicit calculation, that our 
beta functions have the form 



13 gig r) = 2TgR 
(3,{gR) = Dg%^ 



■Ag\- 
Ergl 



Brgl 
-Fgl, 



Cg% 



(45) 
(46) 



where the constants A, B, C, D, E, and F are for now 
undetermined, and computed below. We can already ask 
what is the amount of universality in these coefficients. 
One way to address it is to allow for a class of changes 
in definitions of the renormalized parameters cfr and gR 
such that the new parameters cr^ and are expressed 
in terms of the old ones as 



9r = GgR + HrgR + Ig\ 
<^R = <^R + Kg\ + . . . , 



(47) 
(48) 



where . . . stands for higher-order terms that do not inter- 
fere with beta functions to third order that we are con- 
sidering. Note that a change of the scale of aR is not per- 
mitted since it occurs in the quadratic part, hence is an 
observable. Although one can always consider a broader 
class of changes, this one is broad enough to account for 
changes in definitions of the small- and large-scale cutoff 
and cutoff functions while keeping the structure of p5)) 
unchanged. In particular, the coefficients in Eq. (|43| and 
(|44l) contain some dependence on the details of the cutoff 
function, through the values of the integrals ai and bi. 

One easily finds that the beta functions for the new 
variables is the same as for the old ones (H5|) with the 
change A A/G, C C/G^, D D/G^, and 

B AH + 21 E 2DH + 4:KG 

G2 ' ^G2 G3 ' 

F 2DI-2AGK 

• (49) 

Hence we find that there are the following three invari- 
ant combinations, which we define as: 

AD ^ 4G F + BD-^AE , ^ 
^ = 1^^- (^") 

To see what they mean we consider the value of the 
fixed point I3g{g}f) = of (gS]): 



9R-jr 



2_, 2(2G-A5)^2^^(,3) (51) 



^3 
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a value which is not universal. Then we find that one in- 
variant combination is related to the correction-to-scaling 
exponent by definitioii^: 

^L0^P'^{g},)^-2r + CT^ + O{T^). (52) 

The other two invariant combinations enter into the 
expansion of 

= 2?t2 + (CP - ly +0{t^) (53) 

which will turn out to be related to the amplitude of the 
squared-logarithm, and are universal. 

We now turn to the explicit calculation of the coeffi- 
cients of the beta functions, and of their universal com- 
binations. 

IV. EVALUATION OF BETA FUNCTIONS 

In this section we calculate the coefficients in the beta 
functions. The integrals defined in (gll), (1551) . 

and (j36p are dimensionless numbers of the form: 

Qi — ai{T,'ma), hi — biir^ma), (54) 

i.e., they depend only on r and on the dimensionless ratio 
ma. This is easy to see by the rescaling y — ^ y/{mc) 
which means they can be computed setting m — c — 1 
and replacing — )■ a^rn^f? . Note that for these integrals 
one has —mdm — —ada- 

Hence we will now consider two alternative approaches. 
In the first one (finite ultraviolet-cutoff method) we keep 
finite a in the propagator and evaluate the divergent 
parts of the integrals in an expansion in powers of t. 
In the second approach, close in spirit to the dimensional 
regularization method, we start at fixed r > in the 
glass phase. In that case we find that these integrals are 
ultraviolet convergent, hence one can set a = and com- 
pute these integrals directly in the continuum limit for 
T > 0: 

ftj = ai(T) = ai(T, 0), &i = fej(T) = fei(T, 0). (55) 

The divergent nature of this integrals then implies that 
they admit a Laurent series expansion, i.e., a pole expan- 
sion around t — 0. 

The renormalizability of the theory manifests itself by 
the fact that the coefficients of the above beta functions 
[e.g. P5)) and (011)] will be found to be finite in the limit 
ma — >■ and in the vicinity of the transition r = 0, order 
by order in an expansion in powers of t. 

A. Finite-a method 

For T = ah the integrals ([32]), and 

(155)) contain divergencies when a — )• due to ([T^ . Our 
aim will be to calculate the coefficients close to the tran- 
sition temperature T = Tc (i.e., r = 0) in the form of 



a T-expansion. Inspection of (|43)) and (l44l) shows that 
to obtain the beta functions up to cubic terms we need 
the term ai and 02 evaluated to 0{t) accuracy, while for 
the remaining terms it is sufficient to consider the limit 
T = 0. The detailed procedure for evaluation of integrals 
is given in appendix [Dl Here we only state the important 
results. We mention that a similar method is used for 
the SG model by Amit et al.—. 
The coefficient ai reads 

ai = {2A + tX" + ci [1 + 0{t)] + ©(r^)} , (56) 

where the constant ci is defined in (|D4I) and depends on 
the detailed form of the cutoff function. It is thus non- 
universal. It is important to keep that constant only for 
RG equations beyond gj^. We thus keep it here for our 
two loop calculation and will check that it does not enter 
the final result. We introduced the abbreviation 

X = liiic^m^a^). (57) 

The other coefficient 02 is 

«2 = - ^ - rA - ^ -)- C2 [1 + 0(r)] + 0{r^), (58) 

where the other non-universal constant is C2 defined in 
(ID6|) . 

The evaluation of two-loop integrals 61 and 62 is 
somewhat complicated. The final results are given in 
Eqs. (|D12[) and (jD17[) . The important combination that 
appears in the beta functions now reads 

6i- 862 + 4a? = -A-C3 + 0(r). (59) 

There are two important things to be mentioned. The 
first one is that the A'^ divergence from all summands 
in (j59p vanishes when combined. That is important for 
the renormalizability of the theory and leads to a finite 
beta function. The other point is that the non-universal 
term ciA that appears in 62 and is canceled in the 
combination, leaving the universal coefficient in front of 
gj^ in the first beta function (143]) . 
The last two-loop integral reads 

63 = J [A2 - 2(2c2 + 1)A + C4] , (60) 

with some constant C4. Using mdmX — 2 we finally ob- 
tain the beta functions 

/3gig) = 2Tg-2g^ + ciTg^+g^ (61) 

Mg) = ^5' + (1 + 2c2)Tg2 _ i(ci + 4c2 + 2)g\ (62) 

The beta functions as obtained here are independent of 
the details of the chosen function G{x) to lowest one-loop 
order [the first two terms in (|6ip and the first one in (|62p ] 
but not to the next two-loop order. The apparent one- 
loop universality is only due to our fixed choice of the 
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definition of g in terms of the effective action, wliicli is 
an observable. It is spoiled by any change of scale in g 
[coefficient G in (jTZj)]. or any other change in definition 
of g. 

One thus finds A = 2, C = 1, and D = 1/2 and the 
universal combinations defined in (1501) are 



(63) 



2? = 1/2, C = l, I = F + ^B-E = Q. 



B. Dimensional method 

In this subsection we calculate the integrals by the di- 
mensional method. As discussed above, A divergencies 
obtained using the finite-a method now become poles in 
T. In addition all integrals in the present method become 
TO independent, so only the terms that do not involve 
mdm in beta functions and (|44p survive. Hence, 
from (|43)) and ([44]) we see that the constants in the coef- 
ficients of the beta functions (|45l) and (|46l) are determined 
from the integrals and hi as follows. 

The pole and finite part of the one-loop integrals de- 
termine the four coefficients: 



4ai = -+B + 0{t), 

T 

2a2 = — + E + 0[t). 



(64) 
(65) 



The details of the calculation are presented in appendix 
|F1 with the result 



A 



B 



D = l/2, E 



-4c'i, 
= -2c'^ 



2j 



(66) 
(67) 



where c'l and c'2 are again two non-universal constants 
defined in ^ and (IF^ . 

The two-loop integrals come in two combinations 
which determine C and F as: 

6i-8&2 + 4a? - i- + 0(r°), (68) 

Z T 

6&3 - 8aia2 = - ©(r"). (69) 

T 

From appcndixiFlboth 1 /t^ poles and some non- universal 
1 /r terms cancel in the first combination, which is found 
to have the form (1551) with 



C= 1 

The third two-loop integral is 
1 c' 



6t2 



0(r«) 



(70) 



(71) 



Forming the second combination (|69p one finds a cancel- 
lation of the l/r^ pole and 



F = 2(c'i - c'2). 



(72) 



In summary, one finds the following beta functions 

Pg{g) = 2rg - 2g^ + 4c[g\ + g\ (73) 
Pa{g) - ^5' - 245V + 2(4 - c'^)g\ (74) 

Note that the coefficients A — 2 and _D = 1/2 needed 
for the lowest one-loop-order calculation are the same 
as obtained above in (|6l]) by a different scheme. Again, 
this is because the definition of the coupling constant 
gn is the same in both cases and fixed by the effective 
action. Although the other coefficients, needed for the 
next-order calculation (two-loop) are not the same, one 
can check that the universal invariants 2?, C, and X yield 
the same values as in (|63l) . The two calculations are thus 
consistent. 

For sake of completeness we have performed two ad- 
ditional calculations, both using dimensional regulariza- 
tion. In appendix|H]we have independently confirmed the 
calculation of the present paragraph (and of appendix |F]) 
for a specific choice of the cutoff function G{x) which al- 
lows explicit calculation of all coefficients. In section Ivnl 
and appendices [T] and [j] we have done a different calcu- 
lation using the operator product expansion and a dif- 
ferent cutoff scheme. The beta functions obtained there 
lead again to A = 2, D = 1/2, and C = 1 and now 
B — E — F = which leads again to the same values as 
in (|551) for the three invariants. All four calculations are 
thus consistent. 



V. RENORMALIZATION GROUP EQUATIONS 
AND CORRELATION FUNCTION: RESULTS 

Let us summarize what has been achieved. First we 
have established the following RG equations, in terms of 
the scale £ = — In to: 



0, 



dr 

de ~ 

di 

d<JR. ^ 2 
-,9r 



"^TgR ~ 2gj^ - Brgj^ 

_ 1 

de ~ T 



9w 



Ergl 



Fg%, 



(75) 
(76) 
(77) 



where we recall that t = 1 — T/Tc- Although the more 
general expression are (|45l) and (|46p. from now on, we 
use A = 2, D = 1/2 which, as discussed above, has been 
obtained in several schemes. Here B,E, and F are non- 
universal constant which we found satisfy 



I = F+^B-E = 0. 



(78) 



These equations generalize to next order (two loop) the 
one-loop equations obtained in Refs.— i^^'^'i^. We have 
also clarified their universality to next order in section 
IIIIDl The first equation ([75]) encodes the exact result 
from STS. The first two equations show that the 
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model has a transition at T = Tc- For T > Tc the renor- 
mahzed couphng gn^i) flows to zero, while for T < Tj, it 
flows to a finite value g}^ which continuously depends on 



9r 



0{t% 



(79) 



i.e., the line of nonzero fixed points, shown in Fig. [1] is 
here computed to next order. Its precise value however 
is non-universal. What is universal however is its attrac- 
tive character for T < Tc, together with the value of the 
leading attractive eigenvalue —uj (in the effective-action 
functional space) which defines the correction-to-scaling 
exponent (I52|) : 



-2t- 



(80) 



for r > 0, while it is w = 2t in the high temperature 
phase for t < in the vicinity of Tc- As an example 
of application, we can expect that the dimensionless sus- 
ceptibility fluctuation ratio computed to first order in r 
in Ref.^^ will exhibit a finite size correction as a 
function of system size L. 

We now consider the RG equation for an. From it we 
can obtain the value of the universal amplitude A of the 
squared logarithm, by a simple but non-rigorous argu- 
ment, as was done in^^. Indeed the asymptotic solution 
of mh is 



(81) 



where ctq depends on all details of the initial condition 
and is unimportant as it leads only to a subdominant 
single logarithmic growth. To estimate the off-diagonal 
correlation at a given wave-vector q, one may consider 
the limit of small mass m <^ q and argue that q itself 
sets the scale £* = ln[l/ (aq)] at which one should stop 
the RG. At the same time one replaces a by its effective 
value at that scale, i.e., a <JB,{t}. Hence from (fT5|) one 
writes: 



g^{q) ~ 87r(l 



o87r(l-T)^/3,(g^) 



ln[l/(ga)] 
Hl/iqa)] 



q- 



(82) 



We can now compute the variance of the phase fluctua- 
tions as 



{e{x) - 0{o)y = 2 



(Pq 
(2^ 



{1 ~ cosqx)goiq). (83) 



Using (|82l) and the following estimate of the momentum 
integral: 

d'q , ln[l/(ga)] 
-^(1 - cosqx) 



2n] 
'^-dq 

271 

1 



[1 - Mqx)] 



Hl/{qa)] 



= -^\n^{x/a)+0[\n{x/a)] 

477 



(84) 



we obtain the leading squared-logarithmic behavior 

{9{x) - 6'(0))2 = Aln{x/a) + C'[ln(2;/a)]. (85) 

The amplitude in the above equation is 
A^A{l-Tf(3Ag*j,) 



=4(1 - r)V2[2? + {CV - I)t + 0(t^)] 



(86) 



and using the values of the invariants computed above 
(|63p we obtain our main result 



= 2r2 - 2t3 + ©(t" 



(87) 



In the following section we present a calculation of the 
correlation function from first principles which confirms 
and complement the above more qualitative argument. 

Before we do so, two comments are in order. The cor- 
relation function {[9{x) — 9{Q)f') in the leading order has 
the same behavior as the one in (f85|) . Their difference is 
the thermal correlation function [c.f. that is a loga- 
rithm for the present model due to STS. Another com- 
ment is about the region above the critical temperature 
T > Tc- There ani^) saturates to a finite value at large 
scales. This leads to simple logarithmic growth of the off- 
diagonal disorder averaged phase fiuctuations (|83|) with 
a non-universal prefactor. 



VI. EXPLICIT CALCULATION OF THE 
TWO-POINT FUNCTION 

As is well known, the correlation function can be ob- 
tained from the inverse of the quadratic part of the ef- 
fective action, i.e., in Fourier space 



(88) 



The inversion is in replica space and Tapiq) is defined 
by expanding the effective action (|23p into powers of the 
fields to quadratic order, i.e.. 



r = 



2 ^ 

a/3 



_fq_ 

(2^)2 



To^p{q)9o.{q)dA^q) + 0{e^) (89) 



up to a constant (which encodes for the fiuctuations of 
the free energy) . It is convenient to decompose the replica 
matrix as 



Because of the STS we expect that 



rc(g) 



1 



1 



G(g) 47r(l - r) 



(90) 



(91) 



and one can indeed check explicitly on the expressions 
([24| . ([25]) and ([ST)) that there are no corrections to any 
order in g to Tc{q). This leads to the exact formula 
for the disorder average of the thermal correlation. 
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We now turn to the ofF-diagonal part r(g) which gives 
the second correlation defined in in Fourier space 



Qo{q) = -G{qfT{q) = -{1 - rf 



{Any 



{q^ + m^)' 



rF(g) (92) 



from replica matrix inversion of (|90p in the replica limit 
n — > 042. We now examine its perturbative expansion 
with respect to the anharmonic perturbation ([5]), i.e., as 
an expansion in g, which can be written as a sum 



F(g) = ^r(,)(g), 



(93) 



where F(i) (g) is coming from the corresponding term Fj 
in (ESI). 

The lowest-order term in T{q) is the inverse propagator 



r(o)('7) 



2ti 



(94) 



while the next-order term is trivially obtained from (p4] 
and reads 



(95) 



and it is momentum independent. The first nontrivial 
term is obtained from (1251) and reads 



r(2)(g)--4-.^e-2G(0) 



d2x|e*«"^[2sinhG(a 



(27r)2a4 

- sinh 2G(.x)] +3-4 cosh G(a;) + cosh 2G(a;) | . (96) 

Note that expanding this expression at small q one finds 
that the coefficient of q^ does not yield exactly the renor- 
malized —cr j^/(27r) as defined for convenience in (I57|) . 
The difference however is related to corrections to irrel- 
evant terms in the effective action, and vanishes in the 
limit of zero mass. Indeed the expansion of the second- 
order term of the effective action ([^5)) which produces 
((30)) is somewhat different from the expansion (|96|. The 
reason for that is that the second cosine term multiplied 
by A{x — y,l) from (|^ produced the cosine in (1501) . 
while the first cosine multiplied hy A{x — y,2) from (EH) 
produced the gradient in (l30l) . On the other hand, all 
terms of (|25p contribute to ([M]) . However, in the limit 
771 —> only the term multiplied by A(x — y,2) gives a 
finite contribution that is evaluated below in (j99l) . and 
in that sense the terms ((99)) and the gradient term of 
(|50)) are similar. The finite mass effects are studied in 
appendix [G] 

The previous expressions can be evaluated in presence 
of a small-scale cutoff a in the limit 777 — 0. Using that 
g-G(o) ^ (r77ca)2(i"^) one finds from (ESj) that 



r(i)(g) = o. 

To evaluate (|96p in the limit m - 
note that there is a factor e~'^'^^^^ < 



(97) 

• (and T small) we 
jjj4(i-t) front and 



the integrand can be split in a sum of terms of the form 
e*"-^ — 1. Forp = —2, —1 one can rescale x — >■ xlm and the 
corresponding integrals are convergent in the limit 777 — > 
0, hence the original integrals are bounded by 0(1/777^). 
For "p — \ (and t = 0) the same holds up to a factor 
In 777a. Hence we find that all terms except e^'^*^^-' — 1 
vanish as m — > 0. The limit can be computed by using 
P9p . equivalently written as 



(98) 



The powers of 777 (produced by of e'-^'^^) exactly match 
hence we are left with: 



r(2)(9) =2 



9 



(27r)2a4 
2 



9 



^2 



d^a; 



1 



(1 + x2/a2)2(^-^) 

4"(qa)i-2r^_^^2r(ga) 



9 



- 2t 

qaKx{(ia)\ 



r(2 - 2r) 
r 



aq^^{caq)-\\K^{aq))+0{g 



if 0(09) 
(99) 



Expanding at small aq we find 



r(2)('7) 



-^q'\^{aq)l^--+T{2{:.,,c-\) 
+ 0{g^T\g^q% 



In(ag)] 
(100) 



The cubic term of the two-point function is obtained 
from ()27p . In the limit m — !■ the only term that 
survives contains B{x — y^y ~ z, z — x, 1,1,1). This 
is similar to the simplification which occurs for F(2) {q) 
where only relevant terms in the effective action need 
to be considered in the limit 777 — >■ 0. After trans- 
forming the sum over three different replica indices into 
sums over unrestricted indices using Y^'^^p^ fia, P,^) = 

Eo/.^ 7) - 2E„/(a,a,a) - Ea/^t/K",/?) + 
/(a,/3,a) -I- f{a,l3,l3)], the only term that survives in 
the replica limit is 

r(3)(.)=2^e--(°)^/d2.d2,d2. 



B{x-y,y^ z,z~ x, 1, 1, l)e^{x) [Opix) - 0p{y)] + ... 

(101) 

where . . . stands for all other terms that vanish in the 
limits 771 — > 0, 77. — ?> and for terms that are more than 
quadratic in the fields. After using (|T9l) in the limit t — > 
the previous expression becomes 



r(3)('7)=2^^ / d'^dV(a:-2/) 

^ ' al3 



fix - 2/) = 



x9^{x) [ep{x)-0piy)], 

[{x-yf+a^r' 



' [{y - z)'^ + a'^][{z - x)'^ + i 



(102) 
(103) 
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Doing a Fourier transform one finally obtains 



r(3)(g) 



4g3 
(2^ 



l-c'«-^ 1 



(9 



9ix) = 



(104) 
(105) 



TTln [l + a;2 + |(2 + x^) (.t + V^^Ta)] 
x\J + 4 

After evaluation of the previous integral we find 

r(3)('Z) = f^g'ln(ag)[2(lnc- 1) + ln(ag)] +0(<72). 



27r 



(106) 



To obtain this result one method is to perform the angu- 
lar integral in (|104l) . to differentiate twice with respect 
to q and then to use the following property: 



Axj[qx)- 



■ ln(a; + ..) 



1 

-fiO)\n\ + \nq dz/'(z)ln 
^ Jq 



2 + 0(1), (107) 



where the .. means subdominant terms in the large x 
limit, f(z) vanishes at infinity and xq is arbitrary. This 
formula is obtained by rescaling x ^ x/q followed by a 
partial integration and an expansion at small q. It is then 
applied to f{z) = ^[1 - Jn(z)] and gives (ITMl) . 

We can now add formulas ([Ml), (jlOOp and (|106p and 
obtain 



r(g) = - 



^(g)g^ 

27r 



(108) 



with 



a{q) =a[l + 0{g^)] ~ ^ In(ag) - (g' ~ g'r)[\n\aq) 



+ 2(lnc-l)ln(a(7)]+0(<?3r,gV). 



(109) 



Having obtained (|108p we should reexpress the bare cou- 
pling g in terms of the renormalized one gn. One has 
from (HOl) 



g = g{cmaf^ = .9^(1 + tA -f 2aigR) + 0{g%, r'^gn, Tg%). 

(110) 

Since we want T{q) to order r^, considering that g^ = 
0{t), we obtain 

a{q) =<j[l + 0{g%)] - i ln(ag).g^(l + 2rA + iaign) 

- i9R - 9lr)[ln\aq) + 2(lnc - 1) In(ag)] 

+ 0(g|r,g^r2), (111) 

where to this order we can replace 4ai — > — 2A — ci from 
((56)1 . We recaU that A = In(c^TO^a^). Now, the renor- 
malized coupling gn depends a priori on the product am 
and on the bare value g. However we know that it satis- 
fies the flow equation (pT|) as a function of m. Since here 



we work in the limit m — > it has thus reached its fixed 
point g*j^, hence we must set in the above calculation 

9R^9*n = r+^;^T' + Oir') (112) 



obtained from (|6ip. Remarkably, all non- universal con- 
stants cancel and we are left with 

a{q) = a[l + 0{t')] + l[r^ + + 0{r*)] ln[l/(ag)]. 

(113) 

Another remarkable fact is the cancellation of the contri- 
bution from r(3)((7) with the 0{t) part of the r(2)(g) con- 
tribution. Taken together these cancellations are likely 
to be equivalent to the vanishing of the invariant 1 = 
found in the method using ct/j. One can notice that 
(|94p is the subdominant contribution in the small-g limit. 
Stated differently, the presence of the (bare) disorder 
h{x) in the starting model ([T]), that is characterized by 
the bare disorder strength a present in (jll3p . is not im- 
portant for the leading low-energy behavior of the effec- 
tive action. However, the generated disorder h(x) under 
the RG procedure [contained in terms oc In(ag) in pi3l) ] 
determines the behavior of the correlation function. 
We can now compute the amplitude A. From inserting 



^0 = 87r(l - r) 



cr(g)g2 



(g2 



(114) 



into (j83p . in the limit m — 0, we recover exactly the 
same result as (|57|) for the amplitude. 



VII. RG VIA OPERATOR PRODUCT 
EXPANSION 



Operator product expansion as an efficient tool 
to extract the renormalization constants 



The operator product expansion (OPE) is a very ef- 
ficient tool to extract the RG-functions for renormaliz- 
able field theories. The first to construct a general the- 
ory of renormalization were Bogoliubov and Parasiuk^, 
followed by Hepp^^. They introduced what since then 
is called a R-operation, which subtracts the divergences 
from a given Feynman-diagram, and renders all observ- 
ables, as e.g. correlations functions finite. This was 
done by considering each ordering of the distances in the 
Feynman-integral, the since then so-called Hepp-sectors, 
separately. The R-operation can be thought of as an 
OPE, or Taylor-expansion of a diagram for all possible 
ways to contract the points, from which are retained 
as counter-terms only the divergent contributions, re- 
stricted to the sector in which they diverge. Further it 
could be shown that the R-operation can indeed be inter- 
preted as a multiplicative renormalization, i.e., to intro- 
ducing Z-factors. This was most clearly demonstrated 
by Zimmermann^^, who reformulated the R-operation 
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in terms of forests, i.e., mutually disjoint or included 
sets. An equivalent formulation, which in some respects 
is technically more convenient uses nests. It is this for- 
mulation of the proof of perturbative renormalizability, 
introduced by Bcrgcrc and Lam^'^, which finally has been 
generalized by David, Duplantier and Guitter— , and 
Wiese^ to polymerized tethered membranes, and which 
we will use here. It is a very generally applicable tech- 
nique, which correctly treats distributions, and allows for 
local as well as multi-local divergences. We state the gen- 
eral criterion for renormalizabilit'i^: 

A statistical field theory is perturbatively renormaliz- 
able, if 

(i) the theory is renornializable by power-counting, 
(ii) divergences are short-ranged, i.e., no divergences 
appear at finite distances, 

(iii) the dilation operators commute, 

(iv) there exist an operator product expansion, which 
describes these divergences, 

(v) the divergences of the operator product expansion 
do not have an accumulation point at dimension 0. 
Especially, after subtracting them, the integrand 
has to be convergent when the distances are con- 
tracted. 

This leads us to the theorem of renormalizability 
(i) The renormalized integral 

/ RI{xi, . . . ,xn) 

is UV-finite at r = 0. 
(ii) The renormalized integral, which contributes to the 
connected expectation value of the observable O at 
n-th order, 

O^"^ (zi , . . . , z„) := / R 1 , . . . , f , 

is UV-finite and IR-finite at r = 0. 

(iii) In perturbation theory, the renormalized expecta- 
tion value of an observable is given by 

Or(zi, . . . , Zrn) := V ^-Y- ■ • ' ^rn) ■ 

n=0 

(iv) The subtraction operation R is equivalent to mul- 
tiplicative renormalization, i.e., to introducing Z- 
factors in the standard way. 

In practice, what these two theorems mean for a theory 
like the random-phase sine-Gordon model is summarized 
by the following remarks: 

(i) The microscopic model can be defined without a 
microscopic cutoff, as long as one is below the up- 
per critical dimension (r > 0). It may or may not 
contain a macroscopic cutoff, e.g. the system size, 
(ii) The macroscopic (renormalized theory) is defined 
via perturbation theory. The latter depends on a 



large-scale cutoff L. The choice which is imple- 
mented in the above theorem, is to bound all dis- 
tances which appear in the space-integrals of the 
perturbative expansion by L. If there are strong 
UV divergences, they are to be treated via finite- 
part prescription. 

(iii) Knowledge of the working of the proof of renor- 
malizability, i.e., usage of the operator product ex- 
pansion, is a useful tool to identify and subtract 
subdivergences. The latter can in general easily be 
calculated analytically. Since the resulting (sub- 
tracted) integral is absolutely convergent, it can be 
evaluated in an expansion in r, i.e., for the order 
we need at r = 0. While the OPE is useful to ana- 
lyze subdivergences, it is not necessary to know its 
working to calculate the needed integrals. 

(iv) The theorem then ensures that with this renormal- 
ization all expectation values are finite, for any cut- 
off. 

(v) All critical exponents are universal. 

For the practical calculation, we rely on the techniques 
developed in^"— . The simplest and most pedagogical 
example of a two-loop calculation in this scheme is given 
in^^. An example at four- loop order can be found in^. 



B. Bare Model and renormalization conventions 

We consider the same bare model as in sectionlTTl How- 
ever for simplicity here and below we make the choice 
K = l/(27r), hence Tc — 2 and 

T = 2(1-t). (115) 

The bare action thus is So = H}^^ /T with 

a •'^ "'^ a/3 

where we recall that ^ denotes sums where all replica 
indices take different values. The notation : e*^°*^^^ : de- 
notes normal ordering; implicitly a vertex operator de- 
pending on a single point as e**"'^'' is considered normal- 
ordered, i.e., ^6*^°'"^^) = 1. The correspondence with ^ 
is that go = ^e*^*^"^. The free correlation function for 
.90 = is 

i,Ba.{x)d^{y))=-T8o.0\^\x\. (116) 

Noting the partition function as Z J 'D\Q\Gr^° , the 
effective energy is S = — In Z, defined by expanding the 
partition function to a given order in presence of a back- 
ground field 9, taking the average over the fluctuations 
around this background field, and re-exponentiating, 
keeping only the relevant terms near Tc- It is param- 
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eterized as 



K. 



T 



-E 

a/3 



2T 



27r 

[1 + O(5o) 



E' 

afi 



0[gl) 



T -2t 



0{9l) 



(117) 

(118) 
(119) 



We now compute the renormalized couplings gjf and 
in perturbation theory of the bare couphng g^. We are 
working for r > and in the continuum Umit, taking 



This corrects the quadratic term (in the the Umit of n ^ 
0) as 

SK^p _ 1 



T 



:9q X h, 



1 



(123) 
(124) 



Here and below we are using the prescription that all 
distances are bounded by L, see remark (ii) in section 
I VII Al above. It thus produces the contribution to a^: 



;9o X h- 



(125) 



2. Second diagram (one loop) 



C. One-loop diagrams 

1. First diagram (one loop) 

We use the graphical notation 

„,___0^=e^[«°(^)-«^(^)l . 



(120) 



Further, an ellipse will enclose same-replica terms. Here 
and below we denote SSi = 6si the contributions to 
the effective energy at 1 loop {i — 1,2) and 2 loop {i = 
3,...,6). 

The first contribution is 

— OSl[X) = 



Ci_ / go 
2! V27r 



al3 



.pi[e^{x)-ep(x)-e^iy)+ep(y)] . 



xe 



-2T\n\x-y\ 



(121) 



where Ci — 1 and here and below we use ellipse to show 
which same-replica terms are contracted. The combina- 
torial factor is as follows: a factor of 1/2 from the expan- 
sion of e"''^; the factor Ci = 1 follows from the fact that 
the second pair of replica sums has exactly one choice for 
this diagram. 

This term contains a strongly divergent contribution to 
the free energy (which we do not need) and the important 
sub-dominant term 



Ssi 



1 /5oy 
4 V27r/ 



X -.[{x-y)- Ve^{x) -{x-y)- Vep{x)] 
4 V27r. 



ap 



\x - V\ 



2-2T 



X ■.[veo.ix) -vefi{x)Y : 



(122) 



5S2 = 



■o 



a/37 ^ 



xe 



-T\n\x-y\ 



i[ep{x)+9(<{y)] . 



(126) 



The combinatorial factor is C2 = 2 from the two choices 
to contract the ends. 

Projecting onto the interaction yields 

_ Ss2 ^ x{n-2)Y, 'e'[''°(-)-''-(-)l h , (127) 



a7 



where the basic integral is 

h^:^ [d'y\yr^e{\y\<L)^^ 
Ztt J It 

Therefore 

5(2)5 = (n - 2)gl X h ■ 
where n — ?► has to be taken at the end. 



(128) 
(129) 



D. Two- loop diagrams 

Expanding the partition sum to next order one finds 
the following four 2-loop diagrams. 



1. Third diagram (two loops) 
The first diagram at 2-loop order is 




(130) 



ASc{x)-e^(y)+0^(z)] . 



ap ■''^y^ 

^^-2T[\n\x-y\+\n\y-z\-\u\x-z\] 
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The combinatorial factor is C3 = 3 from the pattern with 3. Fifth diagram (two loops) 

one dot within an elhpse differently colored from the rest. 
This yields 



3! 



The non-trivial integral I3 is, setting y — > 0: 



(131) 
(132) 



\x — z\ 
\x\\z\ 



2T 



e(|x-z|,|a:Uz| <L). 



(133) 

We remind that finite-part prescription is used to define 
this integral. In appendix II II we show that 



h = ^^^^ - 2(1 - + 0(r") . (134) 



= %^^hY^ cosiQ^-Qf,) (139) 
S^'^g = %lh (140) 



The combinatorial factor is C5 = 6(71 — 2)(?i — 3) = 3 x 
2 X (n — 2) X (n — 3) with a factor 3 for choosing the 
middle vertex; 2 for choosing the second white dot in the 
left ellipse - then all vertices are placed; (n — 2) x (n — 3) 
for the replica sums within the ellipses. The integral is 



(2^)2 7x 



6(1x1, |z|, Ix - z| <L) 



WW 



(141) 



2. Fourth diagram (two loops) 



The fourth diagram is 



In appendix II 31 we show that 



T 4t 



(142) 



SS4 



0-- 




(135) 
(136) 



4. Sixth diagram (two loops) 

Ce = 2 is the combinatorial factor of this diagram 
(corresponding to the 2 choices for ordering the 3 vertices, 
up to cyclic permutations). 



The combinatorial factor isC4 = 12(n — 2) = 3x2x2x 
{n — 2) with a factor 3 for choosing the leftmost vertex; 
2 for choosing which of its ends to put inside the left 
ellipse; 2 for choosing the second black dot within this 
same ellipse from the two possible interactions; then all 
combinatorial factors are fixed apart from a factor of n— 2 
for the replica sum in the second ellipse. Setting y — >■ 0, 
we find 



(27r) 



\x — z\ 



e(|a:-z|,|x|,|z| <L). 



In appendix II 21 we show that 



(137) 



2t4t 



16t2 



L^"- (1 - r)L4^ 



O(t0 



(138) 



3! (27r)- 



a/37 

X :Q-^lSc.ix)~e^{z)]^iieffiz)-eiiiy)]^i[e-,iy)^e^{x)] . 



<xJz \x - y\^\y - z\^\z - 



3! (27r)2 2 



■E' 

a/37 



<xJz \x - y\^\y - z\^\z - x\^ 



X -.{(x- z)SIBa + (2 - y)V6'^ + (2/ - x)SIB^ 



2 . 

(143) 



In the analysis, let us distinguish between the resulting 
1-replica and 2-replica terms. We start with the 1-replica 
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term: 



y 



1 rep 



gl ge -3 
3! (27r)2 2 



E' 

Ck./3'y 



'xjz \x - y\'^\y - z\'^\z - xl"^ 



(x - 



3! (2^)2 2 ^ yjj^-2;|T|j^_^|T|^_^|T 



In the last step we did the angular integral (average). We 
need to calculate the following integral: 



ha ~ 



1 



(2^)2 JJ^ \^-y\T\y_,\T\,_^\T 



The calculation in appendix II 41 gives 



(144) 



(145) 



glhaSa^ . (146) 



This yields the correction to Ka/B 

SI<S^ _ C,{n~l){n-2) ~3 
T 3! 2 

We now consider the term off-diagonal in replica-space: 



/ 



2 rep 



3! 



, 1 



3C6(n-2)^' - -.ve^ve^ 



1 



{x- z)- jy- z) 

{2nY JJ^ \:,^y\T\y_,\T\,_^\T 

The integral to be calculated is 

1 f f {x- z)-{y- z) 



(27r)2 JJ^ |^_y|T|y_^|T|^_^|r 



In appendix II 51 we show that 

L6 



6b 



12r2 



+ 0(r") 



(147) 



(148) 



(149) 



This yields the correction for the ofF-diagonal term 

- ^ = ^^^3,3/,,(i - 5^,). (150) 



Taking together diagonal and off-diagonal terms yields in 
the limit of n — )■ 



5K 



(6q) 



SK 



(6b) 



T 



2 ^" 



■^haSaP + /6b(l — <5a,3) 



r 6r 



This gives a correction of cta, 



„3 r6r 

a. = -|^+0(r°) 



(151) 



(152) 



It is an important consistency check that only a gets 
renormalized, but not the diagonal term (temperature), 
as necessary due to the statistical tilt symmetry of the 
problem. 



E. Beta functions to two-loop order 

Summing all one- and two-loop contributions calcu- 
lated above, we find 



9rL' 



-2t 



90 



90 ~9o — 



.90 — 



9 + 6^'^g 
5o 



4t 

r 4r 7" 6r 



O(5o) (153) 
.. (154) 



Here the beta functions are defined as the variation with 
respect to the large-scale cutoff L, keeping fixed the bare 
coupling go- The result, reexpressed in terms of gn, is 



d 
d 



90 



2Tgn^2gl + gl + 0{g%), (155) 
1 



-gi, + 0{g%). 



(156) 



Comparing now with the general expression (|45p and (|46p 
we have obtained the coefficients: 



A = 2, B = 0, C = l, 
D^l/2, E^Q, F = 0, 



(157) 



and we can now repeat the analysis of section IIIII Re- 
markably, the universal invariants in the beta functions 
([50)1 assume, using (|157p the same values (155)) as found 
above. The various methods, while quite different, are 
thus mutually consistent and predict the same amplitude 
A as in (|87p and finite-size-correction exponent uj as in 

(EOl). 



VIII. CONCLUSIONS 

In this paper we have reexamined the random-phase 
sine-Gordon model. We performed a perturbative RG 
calculation in the vicinity of the glass transition tem- 
perature, in a systematic expansion in t = {T — Tc)/T, 
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to the next order (two loop) than was considered previ- 
ously (one loop). We used several different RG schemes 
which yield consistent results. We have obtained the 
scaling equations of the model, i.e., the beta functions, 
given by Eqs. (|75|) . (|76)) . and (|77|) to next order in 
the non-linearity. We elucidated the structure of these 
RG equations which contain several non-universal con- 
stants, and two temperature-dependent universal invari- 
ants. The first invariant yields the correction-to-scaling 
exponent ([5^ which control, e.g. the finite-size depen- 
dence of the susceptibility fluctuations in the glass phase. 
We further calculated the correlation function in the low- 
temperature phase, and found that it has the super-rough 
squared-logarithm form given by Eq. ([85]). Its amplitude, 
which is related to the second invariant in the beta func- 
tions, was obtained to be ^ = 2t'^ - 2t'^ ^ 0{t^). To 
©(t^) it agrees with the one- loop result first correctly 
obtained in^^. The next order Oir^) obtained here is in 
discrepancy with the prediction A — 2t^(1 — r)^ obtained 
in Refi^ by a simple translation to the RPSG model of 
the exact results o:^ based on the fermionic version of 
the model ([I}. The fact that this latter prediction could 
not be correct for all < r < 1, on physical grounds 
and inconsistency with the zero temperature numerics 
was pointed out in Refi^. Here however we find that the 
disagreement occurs already at two-loop level. Hence it 
would be important to perform the RG calculation for 
the model ([T]) directly in a fermionic language and locate 
the origin of the discrepancies between our results and 
the ones from Ref^^S . This goes beyond the scope of this 
paper and is left for future work. 

Another apparent discrepancy can be noted, since a 
calculation taking into account corrections of orders g 
and (which both lead to zero contribution) leads 
to the following result for the Edwards-Anderson order 
parameter—, 



^ei[e(x)-e(o)]^^e-i[e(a:)-e(o)]^ ^ (a;/a)"''^^" 



(158) 



which is different from the one quoted in-® which has a 
temperature independent universal decay exponent equal 
to —4. More studies are then called for to clarify the full 
connection between the model considered in the present 
paper and the one from^S. Note that the predictions 
made in this paper have been compared to a numerical 
simulation, with excellent agreement. The numerical re- 
sults and the comparison is presented in Refi^. 
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Appendix A: Effective action 

In this appendix we derive an analytic formula of the 
effective action up to third order of perturbation the- 
ory. We have not found in the literature a systematic 
derivation of such an expression apart from the result in 
second order perturbation theory in Ref.^**. An advan- 
tage with respect to the traditional diagrammatic per- 
turbation theory is that it gives all the terms in a certain 
order in the coupling constant without need to for deter- 
mination of multiplicative prefactors for diagrams. Our 
aim is to calculate the effective action for a theory de- 
fined by the reduced action (or Hamiltonian in our case) 
S(lp) = So{ip) + gV{(p), perturbatively in g where Sq 
is the quadratic part of the action and V some pertur- 
bation. Denoting by W{J) the generator of connected 
correlations^ 



(Al) 



The effective action is defined as^^ ^{f) = — W{J). 
Using 



J(x) 



5T 



(A2) 



after translating the field ip ~¥ (fi + Xi Eq. (jAl[) eventually 
becomes 



Vx exp 



-S{lp + x) + / <^xx{x) 



5T 



Stp{x) 



(A3) 



Further we introduce T{ip) = r((y9) — S'o(<y5) + InZg, and 
Eq. (jA3l) is transformed into 



1 



Vx exp 



dxxix) 



Soix) - 9V{(p + x) 



ST 



Sip{x) 



exp 



dxxix) 



ST 

Sip{x) 



-gviip + x) 



(A4) 



When deriving Eq. (IA4p we have used that Sq is a 
quadratic action, so it satisfies 



Saif + x) ^Sa{(p) + dx 



SSq 
Sip{x) 



X{x) 



--Sq{(p) 



dx 



SSq 

Sip{x) 



x{x) + So{x). (A5) 



Here one should notice that is important to define F as 
difference between F and Sq in order to avoid further 
complications when solving the implicit equation for F. 
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The next step consists in solving Eq. (IA4[) and extract- and use the cumulant expansion 
ing r out of it. This can be done by iterations. We write ^ 
r in the form i„ _ 1 



n 
ji=i 



where the first few connected cumulants read^ 

{A), = {A) , (AS) 
{A')^^{A')~{A)\ (A9) 
(A^)^ ^ (A^) + 2 {Af ~ 3 (A) (A^) . (AlO) 

Then we easily obtain 

oo 

rM = E3"^«(^) Vii^)^{vi^ + x)). (All) 



n=l 



Next, we get 



(A^)c = 5'(T^^(¥' + X))c + / dxdyG(x - y)^ 



<5r <5(F(^ + X)) 

25- 



Sifiiy) Sip{y) 



-9^{V\^ + x)).-9^ U^<iyG[.-y)'-^^'-^^ (A12) 
7 dip{x) dip[y) 

where we used the notation G{x ~ y) — G{y — x) ^ {x{^)x{y)) ^^'^ the identity 

{V{^ + x)x{x)) = J dyG{x - y)ffl^±^. (A13) 

Therefore, 

In a similar way we obtain 
where we used 

(^(■^ + X)x{x)x{y)) = {V{ip + x))G{x -y) + J dzdtCix - z)G{y - t)^l&±^. (A16) 
Our final expression for the effective action up to third order in the potential V reads 

r(^) .^o(^) - mzo + ,(y(^ + x)> - '^{v^^ + x)). + ^ / dxd,G(. - ,)-^<^(^ + ^^^^^ + 



2 2 7 " ' 5ip{x) 5ip{y) 

_^9'n,3t ^ ^^ 9^ [ A A ^ 5{y{^ + x))5{V\^ + x))c 
+ ^(^ + X))c - y y dxd,G(. - y)^^^ 

[a A A A,n( \rf ,^ S'{Viv + x))S{V{v + x))HVi'P + x)) , 4^ 
+ - y d.d,dzdiG(x - z)Giy - t)^^^-^^ + 0{g ). (A17) 

This formula straightforwardly extends to the case where the fields carry indices, such as replica indices. In that case 
the propagator G carries a double index. These indices can be restored unambiguously using the spatial coordinate 
of the field by matching the field indices with the propagator ones. 

The obtained formula (IA17P agrees with the one of Ref.— to 0{g^) terms. (The formula of Ref.— is written only 
to 0(5^)). We have tested (jA17|) on the sine-Gordon model to two- loop order (i.e., including terms) and found 
agreement with the effective action from Amit et al.—. The final expression (5.1) of Amit et alA contains a typo: the 
argument of the last term in the third line of (5.1) should be y instead of x. 
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Appendix B: Transformation of sums 

When one calculates cumulants of ([8|) one needs to decompose the two sums over pairs of unequal indices, Eqs. (jBip 
and (|B2p . into sums where all addends are sums over all unequal indices [denoted by superscript , i.e., J^a^fi 9i'^j P) — 
5(0;, /?)]. These relations read: 

51 9i.ot, P, 7, ^) = ^ [^(q;, /3, a, 13) + /3, /3, a)] + ^ [^(a, P, a, 7) + g(a, (3, (3, 7) + g{a, P, 7, a) + .g(a, /3, 7, /3)] 
+ ^ 5(a,/3,7,^) (Bl) 

and 

^ g(Q;,/3,7,(5,/i,i^) = ^ 52 + ^ 33 + X! 54 + X! 55 + X! ^6 (B2) 

"7^^ a/3 a/37 a/37(5 a/37(5/i (xfi^S^y 

where 32 = g{a, I3,a, (5, a, I3) + g(a, P,a, f3, f3,a)+ g{a, f3, (3, a, a, P) +g{a, {3, {3, a, f3,a) and 35 = g{a, (3,J,S, The 
remaining terms from (|B2p have somewhat lengthy form and we do not give them explicitly. 



Appendix C: The remaining term of 

In this appendix we quote the remaining term of Fa , that is not important for purposes of renormalization: 

B{x-y,y- z,z- x, ~2, -2, -2) ^ ' cos[6l„(x) + Oaiy) + Oa{z) ~ 0p{x) - 6l^(tj) - O^iz)] 



a/3 

+ 6B{x -y,y~^z,z-x, -2, -1, -1) ^ ' cos[6l„(a;) + 0o,{y) + 6i„(z) - 6p{x) - 0p{y) - 6*^(2)] 
+ 6B{x - y,y - z,z - x, -2, 1, 1) ^ cos[6'Q(a::) + Oaiy) ~ 0a{z) - 0fjix) - 0^{is) + 0^{z)\ 
+ 6B{x - y, 2/ - z, z - a;, -1, -1, 1) ^ ' cos[0a,{x) + 0^{y) - 0p{x) + 0ii{z) - 0^{y) - 0^{z)] 
+ 6B{x -y,y- z,z-x, -1, -1,0) ^ cos[6'a(a;) + 0a{y) - 0fs{y) - 0p{z) + 0j{z) - 0s{x)] 
+ 12B{x -y,y-z,z-x, -1, 1,0)^' cos[0a{x) + 0a{y) - 0p{y) + 0f}{z) ~ 0^{z) - 0s{x)] 
+ 6B{x -y,y- z,z-x, 1,1, -1) ^ ' cos[6la(x) - 0a{y) + 0a{z) + 0fi{y) - 0j{z) - 0s{x)] 

+ 2B{x -y,y-z,z-x, -1, -1, -1) ^ ' cos[6l„(2:) + 0o.{y) + 0c.{z) - 0p{y) - 0^{z) - 0s{x)] \. 



(CI) 



r 



Appendix D: Evaluation of integrals: finite-a method (|43p and ()44p . One should have in mind that we need 

the divergent contributions (when a — ?> 0) from these in- 

1. One-loop integrals tegrals in a power law expansion with respect to the small 

parameter r. The divergent contribution to ai comes 

T .1- J- -n 1 i j-u 1 -J. from the region of integration when the argument of 0(0;) 

In this appendix we will evaluate the unknown mte- . . , . __, .*= , .° , \ ^ 

1,1, • • r ii u J- r m Aix,p] t/Q. (1261) IS around zero. After shifting the 

grals that appear m expressions for the beta functions ^fJi m — cj s 
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variable of integration and using the small-x expansion 
of G{x) given by Eq. we obtain 



ai 



27r 



d^yA{y,l)= / dy 



V 



X {l + T\n[c^m^{y^ +a^)\ + ©(t^)} + f.t. 
- i [2A + rA2 + 0(T2)+f.t.] , (Dl) 



where A is defined in Eq. (j57p . The introduced parameter 
A satisfies 



a < A < (cm) 



(D2) 



and will further serve us to split the divergent part of 
G{x) from the non-divergent one when a — > in integrals. 
In the region (jD2p one can always use the expansion of 
the propagator (|19p . The abbreviation f.t. stands for 
"finite terms" and denotes all terms that do not diverge 
in the limit a — > 0. 

Using a similar reasoning as above, for the other con- 
tribution 02 we get 



a2 



27r 



A'yy'A{y,2) 



dy 







X {l + 2Tln [c^m'^iy^ +a'^)] -f ©(r^)} + f.t. 
1 , ,1 



- |A - rA - ^tA^ + ©(t^) + f.t. 



(D3) 



One may notice that divergent contributions come only 
from the term e*^^^-* of A{y, 1) in (|D1I) and only due to 
q2G(x) f]-QjQ A{y, 1) in (jD3p . The remaining terms in 
A{y, 2) and 1) determine the finite part of the in- 
tegrals. 

For our purpose of calculating the renormalization of 
the effective action up to third order in the coupling con- 
stant g it turns out that finite parts in expressions (jPip 
and (jD3p are important, since they contribute in the beta 
functions. It is important to notice that these finite parts, 
denoted by ci and C2 in Eqs. (|3T|) and (|32|) . multiply the 
renormalized parameters or rg^, so for our order of 
accuracy of renormalization (that is third order in g) it 
is sufficient to evaluate them for t = in (|Dip and (jPSp . 
Here we should have in mind that our general strategy 
of double expansion of the effective action (|37p is in g 
and T, which are assumed to be of the same order. That 
explains why we need constants Ci and C2 to order r'. 

The constant ci , which is a finite term in the expression 
(|D1I) in the limit of r = and for a — >■ is determined 
from the expression obtained from (jPip . 



2tt 



G{y) 



G(2/)-l) = --(2A + ci), (D4) 



taken in the limit a — > and with the propagator G{y) = 
2ifo(m^y2 -j^a'^y This leads to 



ci 



lim 

p-i-O 

1.891 



-4 In (cp) - 4:C^ / dtt 
J p 



,2/fo(t) 



1 



l-8c2 
(D5) 



We should mention that the contribution ?>(? in ci comes 
from the term —G{x) that renders A{x, 1) one-particle 
irreducible. 

Similarly, the finite part in (|D3I) is determined from 



27r 



i2 2 



,2G(y) 



2G(y)-l =-iA + C2, (D6) 



taken in the limit a — for r = 0. This leads to 



C2=-\+ lim <{ In (cp) + c"* 



16c* 



1.611. 



(D7) 



The contribution — 16c^ in C2 comes from the term 
~2G{x) that makes A{x, 2) one-particle irreducible. 
With that, we arrive at the final expressions for ai and 
02 given in the main text, Eqs. and ([55]) . 



2. Two-loop integrals 

We will now calculate the coefficients that stand in 
front of the operators in ([55]) . To achieve that we use the 
procedure described in appendix |E] for evaluating of the 
divergent parts of double integrals. We emphasize here 
that in all terms in we set r = 0, since our purpose 
is to obtain the renormalized action to third order in 
the small parameters g and t, and Fa already contains a 
prefactor g^ . We emphasize here that for simplicity we 
calculate separately the part of B without Bi (terms with 
superscript ') and afterwards we evaluate Bi (terms with 
superscript "), see Eq. ([28)) . The first term of interest is 



6i = 



4 4 

cm 



(2 



^ J d'xd'yBix 



y,x,y,~2,2,2) = b[~bl 
(D8) 



which has divergent contributions from three regions of 
integration, (a), (b), and (c), see Eq. (jE2p . In the first 
region (a) after angular integration we have 



dx 



A 



dyxy 



4x^y 



2, ,2 



(.t2 



',2\2 



4 In — 

4 a 



9 A 
41n^ — 



' A2 



f.t. (D9) 



Combining the contributions from regions (b) and (c) we 
easily obtain 



/•A POO 

b[,^ = - 2c^m^ / d:rx3e2G(-) / dyy/(2, y) + f.1 

Jo J A 



A A 

6 In 81n — In(cmA) -l-f.t. 

a a 



(DIO) 



The contribution that comes from Bi in (IDSP reads 



{2nf 



d^xG{x) / d^y 



'2G(y) 



G{y) 



1 . 
(Dll) 
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It does not contain divergencies and contributes only to 
finite terms. Therefore the final result reads 

h = &'i„ + b'u, - b'l = 5A + + f.t. (D12) 

where A have been defined in Eq. (|57p . 
The second term of interest is 

4 4/- 
C 7T1> 1 

b2 J d^xd^yB{x, y,x + y, 2, 1, -1) = b'^- b'^ 

(D13) 

As in the previous case we consider three regions of inte- 
gration, (a), (b), and (c), see Eq. (jE2[) . In the first region 
(a) after angular integration we get 



b2a = / da; / dyxy 
Jo Jo 



{X^ + a2)2(y2 ^ 0,2) 



1 . A , 2 A „ / a 



-ln-+ln^- + 0^-)+f.t. (D14) 

Combining the contributions from regions (b) and (c) we 
easily obtain 



dx d2/:ry|-x2e2G(-)/(l,y) 
+ e^(^)+^(^)-e«(-)l+f.t. 

A 



-1 + c^m^ 



A. 



dyy 



In. 



In — In(cmA) + f.t. 
a 



(D15) 



The contribution that comes from Bi in (jD13|) is 



cm 



d^xG{x) I d^y 



■My) 



2c^ In(cma) + f.t. 



(D16) 



Using the definition of ci in Eq. (jD5|) after simple algebra 
we finally obtain 



There we used the following result 



A + f.t. (D17) 



eG(3) _ 



2^- — ~ In(cmA) 
4 



(D18) 



where Ci is defined in (|D5p . 

Finally we calculate the integral 



b.^ 



(2 



J d'xd^yx^B{x, x + y,y, 1, 1, 1) = b', - b^ 

(D19) 



Contrary to the previous three cases where the region 
(d) [see Eq. (jE2|] has not yielded divergent terms, in the 
case of (jD19[l it produces divergent terms. First we will 
evaluate (jD19|) in region (a). We have 



<A 

A/a /-A/a 



^G(x)+G{y)+G(x+y) _ ^G{x) _ ^G{y) _ ^G{x+y) _^ ^ 



2ir 



dip 



x^y 



dx I dv 

7o Jo 27r (x2 + 1) (y2 + 1) (y2 + ^2 + 2xy cos (/? + 1) 

'^"d. r\y (l - 

Jo (y^ + l)Viy^ ~ a;2)2 + 2(a;2 + y^) + 1 \ x^ + I 



f.t. 







f.t. 



A/a 



dy 



y 



"2(7/2 + 1) 



In 



A2 



i-y' 



1^ +y^-2(^-l]y^ 



In 2 } +f.t. 



= ln2 _ + f.t. 

a 



(D20) 



Only the first term on the right-hand-side of the first 
line in Eq. (jD20[) contributes divergencies. After using 
the expansion (jlOp . rescaling the variables x ^ ax,y 
ay and doing the angular integration one ends up with 
two terms given in the third line of (jD20p . Only the 



first term gives a divergent contribution when a — )■ 0. 
By doing one more integration over x we end up with 
only one integration over y, see the fourth line of (ID20P . 
Expanding the obtained integral with respect to the large 
parameter A/a followed by integration over y leads to the 
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final result. 

The sum of contributions in regions (b) , (c) and (d) is 
Jo 



X / dyy 
I A 

A 



2c^m^ In ■ 



a 



.2G(a) _ 



dyy^ 



,2G(2/) 



f.t. (D21) 



The contribution from Bi in (jD19P is 

b'i =2|^ / d'xd'yy' [G{y) + G{x + y)] 
+ f.t. = -320^^ In(cTOa) + f.t. 
At the end we get 



,G{x) 



(D22) 



L ul ^ ul 1.11 ^ \ 2 2C2 + 1 
^3 = 63a + ^hcd - ^3 = 4 ^ 

There we used the following result 



A + f.t. (D23) 



/ Azz^ 


e2G(.) „ 1 


J A 





1 



— In(cmA) 
where C2 is defined in (jD7|) . 



C2 + 2 + 16c'' 



(D24) 



Appendix E: Important integrals 

The coefficient in front of operators in T3 contain inte- 
grals of the common type that will be calculated in this 
appendix. They can be written in the following form 



2T.H2„„Pp"G(x)-/3G(2;+y)+7G(y) 



/ = / d'^xd' yyPe' 



(El) 



where the propagator G{x) — 2Ko{m\/x'^ + a^) is ob- 
tained by setting t — into Eq. ([T7| . The param- 
eters in Eq. (jEll) are assumed to belong to the set 
a,f3,j € {0,1,2}, p £ {0,2} that occur in unevaluated 
expressions in Eqs. (|34)) . (1351) and (l36l) . The (logarithmic) 
divergence of integral (IE1[) arises because of the behavior 
of G{x) at X < a ^ m^^. In order to isolate divergent 
from non-divergent parts we split the range of integra- 
tion in (|Eip by a parameter A which is introduced in 
appendix ID] and satisfies ()D2p . We distinguish four re- 
gions of integration 

(a) |a;|,|y|<A, (b) |x| < A, > A, 

(c) |x|>A,|y|<A, (d) \xl\y\>A, (E2) 



and analyze integral (jEl|) in these regions. 

In region (a) one could use expansion ([T9l) for all cor- 
relation functions in (|E1[) and evaluate the integral. We 



will not do it explicitly here for the most general case. 
Particular cases are calculated in appendix [Pl 

In region (b) we can expand the correlation function 
around y: 



G{\x + y\)^G{\y\) 



V 



^G"{\y\) 



1 



\y\ 



■G'{\v\) 



(E3) 



with h = x'^ + 2x ■ y. Then after expanding the term 
Q-l3G{x+y) j^Qj. small h and doing the angular integration 
one gets 

/> A /"oo 

h =(2^)2 / da; / d2/xyi+fe"«(^)+(^-«'^('y) 

Jo J A 

x[l-x'fi(3,y) + Oix^)], (E4) 
where for convenience we have introduced a function 



fW,y) = ^ 



G"{y) + -G'{y) - pG'{yf 
V 



(E5) 



Using similar manipulations as above in region (c), with 
the difference that we expand G{x + y) around x, after 
exchanging the integration variables we get 



4 =(27r)2 / da; 

X [l^x^f{P,y) + 0{x^)\ . 



dyx^+Pye''^''''^+^"^^'^^^y'^ 



(E6) 



While the divergent terms of integral (|E1[) in regions 
(b) and (c) arise only when one of the variables is around 
zero and the corresponding propagator diverges, in region 
(d) both |a;| and \y\ are large. However, their sum \x + 
y\ could be a small number which may in certain cases 
produce divergencies. Therefore, region (d) may contain 
divergencies when |x -|- y| < A and /3 < 0. Changing 
the variables of integration x + y — > x and after using 
expansion (jE3[) one ends up with 

J\x\,\y\>A 

=(27r)2 / dx / dyxyi+Pe-'''5(-)+("+^)G(^) 







A 



X [l + 0{x^)] +f.t. 



(E7) 



The last expression has divergences for (3 = —1. That 
type of integral appears during evaluation of 63, see 
Eq. (p6)) . The case (3 < —1 could be also analyzed, but 
it is of no interest for us. 

We close this appendix by explicitly evaluating three 
integrals. The first one is 



dyy/(/3,2/)=^ 



1 + - + \TL{c^m'^/\^) 

P 

C'(m2A^aVA2) . 



(E8) 
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It could be done by expanding the propagator G{x) = 
2KQ{m\/x^ + a^) around a = to zeroth order, since 
A ^ a. The remaining terms of that expansion produce 
after integration a result which is at least ^ a^/A^. We 
remind the reader that we have already set r = in all 
terms that come with the overall prefactor g^, i.e., in all 
terms that arise from T^. 

The second and the third one can be easily done by 
using the expansion of the propagator ([T9|) followed by 
simple integrations. They read 



1 



■ln(A/a) + 0(aVA2 



(E9) 



and 



A 



dzz^e^GW ^ ^_ [In (A/a) - 1/2] + ^(aVA^). 

(ElO) 

Appendix F: Dimensional method 



In this appendix we calculate integrals (pij) . ([5^ . 
((35)) . and ((36|) by a dimensional method. The main idea 
is to consider t > (i.e., T < Tc) where one can set the 
short-distance cutoff a to zero in all correlation functions 
that appear in the above-mentioned integrals. The log- 
arithmic divergencies contained in the parameter A [see 
Eq. (1571) ] will become poles with respect to r in final ex- 
pressions. The calculation is straightforward once one is 
acquainted with the techniques and ideas presented in 
appendix |D] 

While some results at intermediate steps in appendix 
Id] are calculated using the Bessel function (|17p for the 
propagator, here we show that this is not necessary be- 
cause the universal parts of the integrals come from the 
short-distance behavior of the propagator, that is univer- 
sal. For distances |a;| ^ (cm)^^, it reads 



G(x) = -(1 - T)ln(c^m^a;^). 



(Fl) 



In addition the limiting behavior G{oo) = is necessary 
for the calculation, which is a quite weak assumption. 
Later we will see that we need one more condition and it 
is limj^_j.oo yG'(y) = 0. 

The first term of interest is (PT|) and it can be evaluated 
by splitting the integration range by a parameter A that 
satisfies A <C (cm)^^ [c.f. (jD2[) ]. The divergence arises 
from the region of integration < cc < A, while the 
remaining region x > A delivers a constant . The final 
result can be written in the form 



ai = ^ - c'l + 0(t) 



Similarly we get 



"2 - ^ - 4 + 0{Th 



(F2) 



(F3) 



where again we have a constant that only depends on 
the precise form of the correlation function G{x). We 
should mention that for the special choice of propagator 
([TT]) the constants c'^ and are well defined, however 
their precise value is immaterial for our purposes. 

Further we compute two-loop integrals by the new 
method (they have been already calculated in appendix 
|D]by another method), and we closely follow the notation 
from that appendix. In the following we are only inter- 
ested in the divergent parts of the expressions in the limit 
T — ^ 0. It is convenient to introduce the abbreviation 



9n{A) 

which after using (jFip yields 
ffi(A) = 
52(A) = 



dxx 



-lgnG(a 



1 


■ 1 


c?inn? 


27 


1 


■ 1 
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(F4) 

(F5) 
(F6) 



valid in the limit cmA — > 0. 

The first term is defined in (p4|) . It has divergent con- 
tributions from three regions of integration, (a), (b), and 
(c), see Eq. (|E2[) for the definition. In the first region all 
correlation functions have the logarithmic form (IF1[) and 
we obtain 



b' 



1 1-2 In(cTOA) 

2^2 T 



+ 0(r"). 



(F7) 



Regions (b) and (c) combine into the form given in the 
first line of Eq. (|D10[) . After evaluation one gets 



(F8) 



where /i(A) = dyyC (y)^ . We also used 

£ dyy/(/3, y)=- iyG'(y)|^ - ^h(A) 



2 4 ^ 



(F9) 



+ C'(r°). (FIO) 



where the assumption limj,_>.oo yG'ijj) = has been used. 
The remaining term (IDllI) does not contain divergencies, 
and we finally obtain 

1 3 41n(cmA) -h/i(A) 
2^2 " 27 ^ 27 

Further we calculate the term given in psp . In region 
(a) we obtain 

"■^'jt^-'-^^^ + Oir'l (FU) 

The first two lines of (jPlsp after simple manipulations 
become 



b' 



2bc 



= -32 (A) / dyy/(l,2/)+gi(A)[ai-gi(A)] 

J A 

(F12) 



51 (A)— / d'xGix) 
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while the divergent part of the contribution (jD16l) reads 

b'^^c'^m'^giiA)^ f (fxG{x). (F13) 
2n J 

Combining the previous expressions one obtains 
, 3 4 + 8c'i - 4 In(cmA) - h{A) 



16r2 



16t 



+ 0(r°). 

(F14) 



The important combination that appears in the beta 
function P5|) now reads 

6i-8&2 + 4a2 = -l + 0(r"). (F15) 

There are several important things to mention about the 
last result. First, all non-universal terms connected with 
A and c[ from (jFlOp and (|F14p have canceled in (|F15I) . 
Also there is no 1 /r^ divergence in the combination. This 
is quite reminiscent to the situation we had in the same 
term evaluated by another method, see Eq. ([5^ . 

The last integral we evaluate is defined in ([55]) . How- 
ever, it is more convenient to rewrite it in an equivalent 
form 



(2 



^ J d^xd^y{x + yfB{x,x + y,y,lAA)- 

(F16) 



In region (a) we obtain 

, 1 In(cmA) 
= 67^ + r 



+ 0(t°). 



(F17) 



In regions (b), (c), and (d) now there are two contribu- 
tions. The terms multiplied by x'^ + y^ from (jF16p give 
a contribution that is two times larger than the result 
stated in the first two lines of (|D2ip . The remaining 
term 2a; • y has divergent contributions only in the region 
(d) that equals (|D2ip multiplied by minus one. Overall 
we obtain the same contribution as given in the first two 
lines of (|D2ip . After simple regrouping we obtain 



v. 



Zbcd 



= 2gi(A) [a, - 52(A)] + 4gi(A)i- J d'yy^G{y). 

(F18) 



The divergent part of (jD22p reads 



&'3'=4^5i(A) /dWG(y), 



and finally we obtain 



or^ r 



(F19) 



(F20) 



Appendix G: Two point function at finite mass and 
one loop in the continuum limit a ^ 



In this appendix we study the one loop result (|96p for 
the two-point function derived in the main text. We focus 
on the continuum limit a 0, keeping m finite, which 
exists for t > 0. To recover this limit we must express 
the result in terms of 5, or equivalently of using (|40p . 
We use G{x) — 2(1 — T)Ko{mx), perform the rescaling 
X ^ xj (mc) and get by adding the first and second orders 
in g contributions and (|96p : 



with 



Ao = 
xdx 



(Gl) 



— cWfffl [1 + 2(ai + (G2) 

TT 



-2G(a;) _ ^^-G{x) _ ^G{x) ^ 3 



(G3) 



Here G{x) = 2(1 — t)Kq{x/c), and we note that in the 
combination 



ai -t- ai 



xdx 



-0.0473276, for r = 



(G4) 



the logarithmic divergencies cancel yielding a finite result 
for Aq at r = 0. The momentum dependent part is: 

A(p) = - ^c^m'^gl J d2a;(e'f^ - 1) [2 sinh G(a;) 

-sinh2G(a;)] (G5) 

Let us now consider its limit for large p — q/m limit 
obtained by rescaling x — y/p. Using 



lim smh[2G{y/p)] = 1/ h/'^-^) 



p—^oo 



lim smh[Giy/p)] = 0, 



(G6) 
(G7) 



one finds 



A{p) 



TT 
1 



-2 2 2 2-4t 

-c m gj^p 



ydy{Jo{y) - l)y-^^'--^ 



IQ 

24t-3t 



-r'^m'^n'^ „2-4r ^ r(2r 1) 

^cmgj,p r(2-2r) 



1 



2 2 2 

cm gii 



P 

16r 



1 



p^ \ii{pc/e) +0{t) 



(G8) 



(G9) 



Let us now comment on the various regimes for T{q) as 
a function of q. By comparing with (|100p we see that 
the In q behavior (that leads to super-rough correlation) 
extends from the region g '-^ 1/a up to the region q ^ m. 
However for very small q, q ^ me^^^^-^'^^ the coefficient 
-—2TTT{q) / q^ saturates to 1/8t (plus a finite part). The 
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pole in r has precisely the value obtained in (|39|) together 
with ([55)) and (p7|) by considering the limit g — ?> first. 
Finally the above result allow to determine the correla- 
tion function Goiq) = /('?/'7i)/m^ as a scaling function 
of q/m using (|92t . 



Appendix H: More on two loop integrals: 
simplifications and explicit evaluation for a simple 
cutoff function 



In this appendix we calculate the integrals ([3T|) . (132 
and ([M| - ((5S|) using a simple cutoff function. 



1. Simplifications 



First we show that the terms linear in G{x) in all inte- 
grals, which make the effective action a sum of one par- 
ticle irreducible graphs, cancel in universal invariants in 
the beta functions. We also provide simpler expressions 
for evaluation of two-loop integrals. 



Let us start with one-loop integrals. We define: 



2tt 



V 



pG{y) ~ 1, 
- pG{y) - 1 



(HI) 
(H2) 



so that oi 



,(1) 



and 02 



,(2) 



Here we denote f 



J d'^y. Analogously we define expressions without the 
term linear in G{x): 



2 2 

2-K 



2tt 



y 



.PG{V) 



(H3) 
(H4) 



This changes only the finite parts in the one-loop inte- 
grals since the integrals linear in G{y) are finite. More 
precisely ([M)) and (|55t are changed into 



with B 



4ai : 
2a2 = 
B 



4a 



(1) 



A 



B + Oir), 



2a. 



(1) 



T 

-+E + 0{t), 



(H5) 
(H6) 

E - 



4%rLG(x) and E 

42g^ L x'^G{x) (in other words i = A and £» = £>). 

For the two-loop integrals we will use an expression 
equivalent to 



B{x,y,z,a,b,c) 



^aG(x) 



^cG(z) _ ^ 



.cG{z) 



-cG{z)-l . (H7) 



Then one can rewrite the needed combinations of two-loop integrals: 



bi - 862 + 



/ x-\-y-\-z—0 



■8 



-2G(z) 



,2G(a:) 



^2Gix) 



.2G(y) 



+ (a(2V+2af^a(-^)-t-4(a(i))2 

finite (H8) 



, .(2).(1) , .(2).(-l) , .(l).(-l) 



the same expression without the hats and the additional finite part being exact [by definition from (IH7p ]. Comparing 
with ([68)) it means C ~ G, i.e., C can be computed discarding linear terms in G(y). The second combination (|69)) 
can be rewritten: 



663 — 8aia2 =6 



(2^)2 



x+y+z—O 



24dia2 — 80102 -|- i 



.G{x) 



27r 



1 



gG(y)_^ e«W-l 

a2G(a;) — aix'^G{x)\ + finite 



(H9) 



I 

again the same expression without the hats, the terms linear in G, and the additional finite part being ex- 
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act. Comparing with ([69)1 one sees that now the coef- 
ficient F is changed, i.e., F = F + AD^^ G{x) - 
2A "]^^ J^x^G(x). However the combination 

F + BD-^AE = F + BD-^AE (HIO) 

does not change and since A and D are also unchanged 
the three universal ratios (l50t are unchanged. Hence one 
can just suppress the terms J^G{x) and J^x^G{x) since 
they cancel in universal quantities. 



2. Calculation with a simple cutoff function 

We consider the continuum (or dimensional) limit 
scheme where a — and r > 0. Until now we have 
used the Bessel cutoff function G(x) = 2(1 — T)Ko{mx), 
and noted that all factors mc in the definitions of the 
integrals can be set to unity by rescaling x — x/lmc). 
This is equivalent to using G{x) = 2{1 — t)Ko{x/c). Here 
we will test further the universality with respect to the 

I 



choice of the cutoff function by considering 
eG(^) ^ e(a; < l)x2^-2 + Q(x > 1) 



(Hll) 



such that G{x) is continuous and has the same logarith- 
mic behavior at small x as the Bessel function cutoff and 
vanishes at infinity. For convenience we also suppress the 
tilde (and set m = c = 1 in all definitions). 
The one-loop integrals are: 



ai 



a2 



dxx{x'^ 



Axx {x 



1 1 



so that 

^ = 2, B 



-2, D = l/2, E=-l/2. 
a. Coejficient F 



(H12) 
(H13) 

(H14) 



From (jH9l) and (p^ the coefficient F can be extracted as 

- + 0{l) = 6 f a;2(eG(") - l)(e^(y) - l)(e<=(") - 1) + 24aia2 - 8a,a2 = - ^ + - 



(H15) 



where we have used (jH12|l and computed the one loop integral 0,10,2 = ^ Jq dxx^{x ^ — 1) = l/(8r), keeping only 
divergent parts. We have defined 



^1 



1 



2{x^ + y"^ + z'^)ix 



2t-2 



2r-2 



1)(Z 



2r-2 



(H16) 



(^TI")^ Jx<l,y<l,z<l 

To compute this integral (which is ultraviolet convergent) we use symmetry to restrict to the domain 0<y<x<z<l 



and write y = wx and < it; < 1 and z 



-y^ + 2xyc with c = cos(0). The measure is then 12a;^wda;dt(;dc/ \/l — 



with c G [—1, 1]. The condition x < z < 1, i.e., x^ < z^ = x'^{l + 2cw + vo^) < 1 implies that 2c + w > 0, hence we 
use that 



. . = 27r 



x<l,y<l,z<l 



dc 



1 VT 



C^JO 



1 pl/Vl+2cw+w'-' 

dwl2w / dxx^Q{w + 2c> 0) . 

Jo 



(H17) 



After integration over x for < X < 1/VT+ 2cw + 
one finds 



with 



1 

2^ 



dw 



' dc 



VT 



^gr{w,c)Q{w + 2c> 0) 



(H18) 



..Kc) = l ^[f 7^ + V o(r°) (H19) 

That would naively give A^i ~ A/t but there is a pole at 
w = 0^. To treat this pole one checks that replacing in 



dHTl gr = gr - <?? with 



16(r - l)2w2^-i 
,?(z.,c).-^-,J^e(0<.<l) (H20) 



one finds a finite integral ~j4/rasT— >-0. Hence we can 
write A^i = TVf -I- Nl'^^ where the contribution of the pole 



27 



Jo 



dc 



1 4(r-l) 



47r - 3 X 4^Bi t 



2' 2/1 



27r 



3t2 (2r2 + 5r + 2) 



(H21) 



Here a = 



3V.<^'(l/3)-V/^'(5/6) 



reg _ 



and the regularized part is 

_J_/"'du;/" dc 1 8c 

^2^ Jo ^ J-^u/2 t2cw + w^ + 1 

+ 0(r°) = -- + O(r0) (H22) 

T 



where we have checked that the remainder ©(t^) is inte- 
grable. In total we have A^i = I/t^ — 9/(2r) which gives 
F = 1/2 and the combination F + BD - ^AE = using 
(|H14p . Hence we find that the universal combination is 
X = 0, in agreement with the other calculations. 



b. Coefficient C 

To treat the other two-loop coefScients one should 
group the terms so that the integral is ultraviolet finite 
for r > 0. This leads to: 

+ O(r0) = bi-8b2 + 4.al^M + Q+\-- 
2t t 



(H23) 



with 



M 



-2G(z) 



.^2G{x) 



,2Gix} 



Q = -8- 



(27r)2 7a:+y+z=0 



^2G{x) 



^2G{y) 



-G{z) 



■My) 



-2G(z) 



1 



,2G{x) 



(H24) 



1 



(H25) 



After symmetrization on the arguments x,y, z one ob- 
tains (with X + y + z — implicit): 

M + Q = Ni+N2, Ni= [ fi{x,y,z), 

J x<l,y<l,z<l 



N2 = ^ I f2{x,y,z), 

' x<l,y<l,z>l 



(H26) 



since when two or more distances are larger than 1 the 
integrand is identically zero. The functions /i and /2 
are easily obtained using Mathematica. For the sector 



X < l,y < 1, z < 1 we perform the same manipulations 
as above. We find again (jHlSp with now 



?t(w,c) 



-8 {w [w^ - 2c''w -I- 6c^w + 3c {w^ + l)] + l} 



Tw {2cw -I- -I- 1) 



0(r°) 



(H27) 



The pole at w = can again be treated via a substraction 
with now: 

gP{w, c) = \ e(0 <w<l). (H28) 

t{t + 1 ) 

To see this we note that g-r can be split in three terms: 
(i) the coefficient of w'^+^t which gives the above (ii) 
the coefficient of w"'^^^'^ which turns out to simplify into 
an expression yielding no pole in w (due to a cancellation 
of poles between M and Q) (iii) an expression regular at 
w = 0. Hence we have again iVi — Nf + N^^^ where the 
contribution of the pole is: 



^ 2(r- 1) 47r - 3 X 4^Si (r 



1 1^ 

2' 2) 



27r 



3r2 



1) 



1 



2-a 



T 

and the regularized part is 



+ 0(r°). 



(H29) 



1 

'2^ 



1 .1 
dw \ 

J-w/2 Vl 



dc 



8 (2c^w - 2c^w + cw'^ + 2w) 



{2cw + w'^ + 1) 



0(r°) 



(H30) 



In total we have iVi = — I/t^ -|- 5/(2t). Next, one easily 
checks that N2 is finite, as the expansion of the inte- 
grand {gr) up to 0{t) at small w is regular. Hence from 
(IH23P and below we find cancellation of l/r^ poles and 
the residue C ~ \. This is again in agreement with the 
other calculations. 



Appendix I: 2-loop integrals for the OPE method 
1. Integral /s 

The first integral at 2-loop order is 



1 



2T 



e(|a;-z|,|a;|,|z| <L) 



(II) 

We remind that this and other similar integrals are 
defined via finite-part prescription. Integral (HTj) is 
a counter-term of the theory, which has itself sub- 
divergences. It can be calculated by brute force, but the 
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lzl> 1 




FIG. 3. The geometry used for J3 



task is simplified by recognizing that the structure of the 
renormalization group dictates these sub-divergences, or 
equivalently lower-order counter-terms. After subtrac- 
tion of these counter-terms, which themselves are easy 
to calculate, the rest can be written as a convergent in- 
tegral and be calculated. This is the road we will follow. 
To show how this works, consider the OPE-coefficicnt 
associated to I3, which is the integral of three interac- 
tions projected onto a single one (not writing the replica 
content), which we note as (not yet setting y = 0) 




■0 = 



\x - y\\z - y\ 



2T 



(12) 



Subdivergences occur for x ^ y or z ^ y. Let us consider 
X ^ y, keeping in mind that a similar term appears for 
z y. Then, we use the OPE to extract the small- a; 
behavior: 




-O) 



■-Oz m- 



■o) 



X (VjVfc6>(2/) 




X (V,9{y)VkO{y) • 



^Ay)^kO{y) 



--Oz 



-o) 



(13) 



The above symbolic notations mean that we first consider 
the OPE of the upper two interactions, and projecting in 
the order of their appearance: on the identity operator 
at position y, ly (termed relevant counter-term below) 
times the contraction of ly with the interaction; then on 



'Vj6{y) (termed relevant odd counter-term below) times 
the contraction of V j9(jj) with the interaction; then on 
VjVkdiy) times the contraction of the latter with an 
interaction; and finally on the interaction itself times the 
projection of the interaction with the final interaction, 
termed marginal counter-term, or RG below, since this is 
the repeated counter-term from 1-loop calculations. The 
coefficients in the order of their appearance are 




ly ^\x-y\ 



-2T 



Ozl •---O ) = 1 




V,0{y) 



--Oz •----o 



VjVkOiy) 



'2i{x-y)j 
\x - yp^ 

^ iT{z - y)j 
\z-yr 

i{x-y)j{x-y)k 



\x-y\ 



2T 



(VjVfc6l(2/) •---O^l •---O 



iT[djkiz - yf - 2{z - y)j{z - y)k\ 




\z-y\' 
^Ay)^kO{y) \ = 



(14) 
(15) 
(16) 

(17) 

(IB) 

(19) 



-2{x - y)j{x - y)k 



\x - y\ 



2T 



{V,0iy)Vk0iy)9---<^z\ #---0) 

^ -T^iy-z),{y^z)k 

\y~z\^ 



(110) 
(111) 



Some remarks are in order: In a first-principle calcula- 
tion, the replica content would have to be written explic- 
itly. Thus, e.g. the factor of 2 in (jl6| comes from the 
2 possibilities to expand either the left our right replica. 
In the effective action calculation presented earlier in this 
article, the first three terms of p3|) are absent, since they 
correspond to 1-particle reducible diagrams, which are 
automatically subtracted. 

This allows us to give the list of counter-terms. Note 
that if we only know that such an expansion exists, but 
do not know the coefficients, Eq. p3|) can be obtained via 
simple Taylor-expansion. (That is actually what we did, 
except to check the working of the OPE on an example.) 

Counter-terms for integral I3: To simplify notations, 
we set y — 0. First, the relevant counter-terms, minus 
the product of p4| times (|I5|) . since we want 73 plus the 
counter-terms to be finite, (plus an analogous term with 
X and z exchanged) are 



1 



i2ny 



1 



1 



\2T 



2T 



e(|z 



.(112) 

This expression is zero due to analytical continuation. 

Second, the odd relevant counter-terms, minus the 
product of (|I6)) times (|I7|) (plus an analogous term with 
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X and z exchanged), 
2T 



TC.odd 



(2n) 



2T 



\z\^\x\ 
xe(|x|,|z| <L) 



1^1 



2T 



(113) 



These terms are zero due to analytical continuation, and 
zero due to parity. The marginal counter-terms reads 



rC,RG 



T 



i2ny 

J2 .,2 



+ 2{T - l){z ■ xY 



■e(|x| < \z\ < L) 



x'^ z"^ + 2{T - l){z ■ xf 



■e(|z| < \x\ < L) 



(114) 

We have written two contributions: the first ^ is 
the repeated counter-term from RG, equal to minus 
(|/10P X (j/ll[) (plus an analogous term with x and z 
exchanged). Note that we have put 8-functions to re- 
strict the counter-term to the sector in which it is needed 
to subtract the divergence. The second contribution is 
(j/Sp X (|/9p (plus an analogous term with x and z ex- 
changed). It is zero due to radial integration, but makes 
the integral absolutely convergent (it can e.g. be put on 
a computer). 



Consider J3 ~ L-^[I^+i.^ + jg -± 1^ 
Using the mapping theorem of appendix [Jl to map onto 

\x — z\ = L = 1, using z = X — 1. we can write with the 
respective counter-terms regrouped in one line 



xc.ici , xc.odd , J 

or ^ ''3 + -'3 _ + ''3 



cRGi 



4r 



= L 



d 



dL 
1 

2^ 



xc.rcl . 7-c, 
'3 + ^3 

I T-c.rcl . 



3dd 



c,RG 



L=l 



rc.odd 



max(l,|a;|,|l-a;|)-^^ 



c,RG 



L=l 



\2T\ 



1-x 



1 



2T 



1 



|2T 



|2T 



|xr- |l-a;|- 
2T[x ■ (1 - x)] max(|a;| 



x\) 



-4t 



|x|2r|l-2;|2 
2T[x- (l-a;)]max(|a;|,|l 



x\)- 



\2T 



-T 



x^{l 



2{T-mi-x)-xY 



\2T 



(1 



xe(|a;| < |l-a;|)|l 



-4r 



x2(l-.T)2+2(r-l)[(l 



xe(|i 



X < X 



(115) 



J- 



normal 



J: 



anomal 



There are normal and anomalous terms, where the latter 
are those for which the difference between the above max- 



functions matters, even though in the limit of r — ^ they 
all become 1. 

The normal term is 



J" 



)rinal 



1 

2^ 



1 



1 



1 



\2T\ 



\2T 



\2T 



|2T 



2T[x-{l-x)] 2T[a; • (1 - a;)] 



\2T\ 



|2T 



T-^(^-)^+?i^r^l[(^--)--]^9(N < II - .1) 



T 



|x|2T(l_x)4 

x^l ~ xY + 2{T ^ ~ x) ■ x]^ 



\l-x\^T^i 
X max(l, |a;|, |1 — x\)~^'^ . 



e(|i-x|<|x|) 

(116) 



Since we have constructed the counter-terms such that 
this integral is convergent in each sub-domain, one can 
take the limit T — > 2, i.e., r — > 0. This yields, up to 
terms of order r, 







(117) 



since in that limit the integrand vanishes identically. 

We now turn to the anomalous terms. Those come 
from the regions where the divergences do not cancel ex- 
actly, either x — >■ 0, or a; — >■ 1. Due to symmetry, we 
consider a; only (this gives a factor of 2). We write 
the different contributions as follows: 



J: 



|„omai ^ (|7T5D_(|7T6D 



janomal_^ja„omal^0(^-) (j^g) 



The first term is 



ranomal 
'3a 



= -2— 
2tt 



1 — max(l, |a:|, |1 — a:| 



|2T 



xe(|.T| < 1/2) , 



(119) 



but one could have taken a smaller number than 1 /2 for 
the cutoff. This gives 



J; 



anomal 
3a 



.2^ 

2tt 



- 1 



•|2T 



e(|f| < 1/2) 



xe(x- 1 < 0), 



(120) 



where we have made explicit that x is a vector, see Fig. 
[3l Going to complex coordinates gives 



J: 



anomal 
3a 

2 

2^ 



2 

2^ 



(l-a;)-2-(l-a;)- 



- 1 



[xxy^ 

X e(|x| < l/2)e(5R(x) < 0) 
2t(x + x)+ At'^xx + 2t(2t -|- l){x'^ + x^)/2 
(xx)"^ 

X e(|a;| < 1/2) e(5R(x) < 0) + . . . 

I /•37r/2 p 



2 / dxa;^^-^ — 
/o 2^ 



Tr/2 



2Ta;(e''P + e" 



2^2 



+ t{2t + l)x''{e^"^ + e-'"^) + 



(121) 
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since the only (logarithmically at r = 0) diverging terms 
are the term ~ At'^x'^ — ?> ^i''') after integration, and the 
last one, which vanishes due to angular integration. 
The interesting term is 



J: 



anomal 
3b 



X [max(|f|, \l-x\)-^^ - max(l, \x\, 



4T f x-{l- x)[\l-x\ 



27r 



-4t 



-1 



•|2T| 



1-f ■ 



X e{\x\ < i/2)e{x- 1 > 0) 

4T r ^{x + x- 2xx)2t{x + x + ...) 

'2^X; ixx)'^{l- X){1- X) 

X e(|2;| < l/2)&{^{x) > 0) 
4Tt f^^^ 



27r 



d(/i / dx x' 



,4t-3 



ir/2 

(xe*"^ + xe-''!' - 2x'^){xe"'' + xe'"!' + ...) 



(1 - xe*'^)(l - xe-^'t') 



ATt ['"I'^ 



2tt 



l-n/2 Jo 

-T + 0{t). 



(122) 



The remaining anomalous terms only diverge logarith- 
mically at small x, thus when expanding the factors of 
max(...)~'^'^, this yields an additional factor of t\x\, en- 
suring convergence; thus these terms can be neglected. 
Therefore, calculating (jI14p analytically, and using (jllSP 
we conclude that 

T rC.R.G I 1 Tanomal 



At 

2 t4t 



'3b 



2(1-t)— + O(t0) . (123) 



At 



2. Integral I4 



h 



{2n) 



\x — z\ 



Q{\x-z\,\x\,\z\<L). 



(124) 

We follow the same strategy as for , first identifying the 
counter-terms. The relevant counter-term to be added is 



J- A 



1 



1 



{2nfJJ, \x\ 



\2T 



e(|a;-z|<L). (125) 



We can again add the following (relevant) odd (i.e., van- 
ishing) counter-term 



3dd 



T 



(27r)2 JJ^ \x\^T,2 



e(|x|,|z| <L). (126) 



This time, there are two marginal counter-terms. The 
marginal counter-term for x — comes from the inser- 
tion of (fmj) : 



rc,RG,l 



T 



i2n) 



(r-2)(x-z)2+x2z2 



xe(|x| < \z\ < L) 



(127) 



The marginal counter-term for z — > comes from the 
sub-divergence (|126p : 



c,RG,2 



1 



1 



(27r)2y,A IxFlzF 

Consider now the combination 
d 



e(|z| < |x| < L) (128) 



J. 



dL 

1 

2^ 



rc.rcl I rc.odd 

,1 1 /I 



tC.RG.1 , rC,RG,2 
+ 



1 



max(|x|, 1, |1 — x\y 



-4r 



|2T 



-T 



X ■ {1 — x) 



\2T\ 



l-x[ 



max(|x|, |1 — x| 



^ (r-2)[x.(l-x)]^+x^(l-x)^ ^ 



2x 



2T\ 



xe(|x| < |i-x|)|i-x|-*" 
^ .e(|i-x|<|x|)N 



x|J |l-x[ 

anomal 



jrnormal _|_ j. 

The normal contribution is 



(129) 



rnormal 



2tt 



|2T 



T 



X ■ {1 — x) 



\2T\ 



1-x' 



^(T-2)[x. (l -x)F+x^(l^x)^ ^^1^1 ^ 1^ _l 



1 



ixrii 

0{t) . 



2|x|2T|i_a;|4 
^e(|l-x| < |x| 



max(|x|,l, |1 - x\)-^'' 
(130) 



The integral is convergent, since the counter-terms where 
constructed in order to cancel exactly all sub-divergences. 
Taking the limit of r — )• 0, one finds that the integrand 
identically vanishes, which shows that the expression is 
0{t). 

There are two potentially anomalous terms, 



ranomal yanomal 

•^4 — •J4a 



J: 



anomal 
46 



(131) 



For the first term, a calculation identical to (|I19p shows 
that 



anomal 



1 

27r \x 
0{t) 



^[max(|x|,l,|l-x|)-4-_l] 

(132) 
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The second term is, up to a prefactor of 1/4 identical to 
(II22]). and reads 



j-anomal ^ f 2^ ' (1 2?) 



Mb 



2^ .L |f|2T|i_f 



e(|x| < l/2)x 



[max(l, |f I, |1 - xl-'^^) - niax(|f|, |1 - f 1)-"^] 

(133) 



Tanomal 
'^3b _ 



4 4 

Calculating explicitly the integrals (|I27I) and (II28P yields 



c.RG.l rC-RG,2 , 1 ranomal 



T 7-C.l\,«^,l rl 

J4 — ^-'4 ^ 



— 



(1 - r)2L4- L4r _ 



4r 



16t2 8r2 8t 



+ ©(t") . (134) 



We remark that the contributions proportional to ~ 
(1 — r)2 from I3 and J4 cancel. We will see later that the 
only terms which appear in the RG-functions come from 
the anomalous terms. 



3. Integral Is 



The integral I5 is 

= J— f f Q(I^U^U^-^I <-^) (135) 



(2^)2 



\x\T\z\T 



Clearly, the marginal counter-terms are subtracted by 
h - II = 

1 f f (d{\x\,\zl\x - z\ < L) - @{\x\,\z\ < L) 



{2n) 



(136) 



Using the mapping prescription, we find 

1 /■ max(|a;|,|l-x|,l)-''^-max(|2;|,|l-a;|)-'*^ 

= 0{r) (137) 

Thus there is again no genuine contribution. This is not 
astonishing for a bubble-chain. The final result is 



/5-/| + 0(r") = ^ll + 0(r°) 



(138) 



It has two subdivergences, due to 1-loop counter-terms 
for the coupling g: for x — and for y — z ^ 0, 



^6a 



-1 



(27r) ,,,, 
1 



(140) 



\x - 2/r \y - z\ 

1 {x - yf 



\z - y\^ \x - y\ 



2T 



Q{\z-y\ < \x^y\<L) 



We have explicitly written the subdivergence (first fac- 
tor) times the remaining term (second factor) times the 
restriction on the sector (third factor, 0-function). Note 
that we have used our freedom to put the second factor 
at a point of our choice. 

We now note that (i) when combining the integrands 
of Iga + -^60^^' there are no subdivergences, and (ii) the 
integrand vanishes in the limit of r — 0. Therefore, 



ha 



^c,RG 
'6a 



finite 



and 



jcRG 
'6a 



O(t0 



^6r 



-O(rO) 



(141) 



(142) 



5. Integral let 

The integral /gb is 

1 f f {x- z) ■ {y- z) 



'6b 



(2^)2 JJ^ \:,-y\T\y_,\T\,_^\T 

xe{\x-yl\y-z\,\z-x\<L). (143) 



It has a sole subdivergence, when x — y ^ 0. It is sub- 
tracted by 



^c,RG 
'6b 



1 



(2^)2 



\x-yr \y-z\ 



{y~z? 

2T 



(144) 



X <d{\x-y\ <\y-z\<L) 



We have decided to use this specific form of the counter- 
term (for later convenience) but a symmetrized version 
would also be possible. We claim that /gb + Iq^^^ — 
©(t"). To prove this, we set y — 0, vary L, and map 
onto |a;| — 1 



4. Integral ha 

The integral /ga is 

1 r f (x - zf e{\x - y\,\y - zl\z - x\ < L) 



\x - y\'^\y - z\^\z - x\^ 



(139) 



" 2^ 



L = l 



z • (z — 1) max(|z|, 1, |1 — z\ 



e(|z| > l)max(|z|,|l-z|)- 

|,|2T-2 



(145) 
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The integral is finite, thus one can go to the critical di- 
mension. This yields 

Ld 



dL 
2^ 



L = l 



z-iz-1) e(|z|>i) 



(146) 



The integral can be split into two parts, Jq = + , 
calling the part where \z\ < 1, and the other term the 
part for \z\ > 1. We get 

ZTT / z^\z — Ir 



i^/d^.e(N<i) 

d^zQ{\z\ < 1) 



1 1 

42^ 



\2z-l\^ - 1 

Z2|z- 1|2 

(2z-l)(2z* -1)-1 



zz*(z — l)(z* — 1) 

where in the last line we have introduced complex coor- 
dinates. This gives 

1 1 

z(l — z*) z*(l — z) 

(148) 

Since \z\ < 1, Taylor expansion can be used around zero. 
It shows that there are only terms which vanish upon 
angular integration. Therefore 

. 



^a< = 4i/d^^e(N<i) 



(149) 



1 



z(l — z*) z*(l — z) zz* 
Using that the Taylor-expansion for z — > oo is 



11 1 1 

z 1 — 1/z z 



1 - z 



1 



1 



(150) 



(151) 



we conclude that the terms of order l/(zz*) cancel, and 
the remaining terms vanish upon angular integration. 
Therefore 

J6> = . (152) 

Putting the pieces together, we conclude that 

hb + /gf = finite . (153) 

We checked that one can also map onto |a; — z| = 1, 
yielding 



Ji 



6b 



Ld 
dL 

2^ 

2^ 



L=l 



tc.RG 
hb 



(154) 



l-(l-x) e(|a;| < 1) 



1 



2xx* 



2(1 - x){l - X*) 

1 e(|a:|<l) 



2xa;*(l - x){l ~ x*) 



xx* 



= 



using again the Taylor-expansion method, this time sep- 
arating into |a;| < 1 and |a;| > 1. Thus, renaming the 
variables, 



6f) 



^/e^f« + O(r0) 



(27r)2 

6t 



i^Se(N<|z|)-fO(r°) 



L 



\2t' 



+ O(t0). 



(155) 



We note a consistency relation between integrals /ga and 
/gfc: Rewriting the numerator of (jI39[) as {x — z)^ = (a; — 
?/)2-|-(?/ — z)2-|-2(a; — y)-(?/ — z), we deduce that /ea = 2/6&. 



This is indeed satisfied by our results (|I42[) and (|I55 



Appendix J: The conformal mapping theorem 

In this appendix we discuss a convenient method to 
extract the divergence of an integral, known as the con- 
formal mapping theore m^^'^^'^^ . This method we use to 
extract the 1/r contribution of the 2-loop integrals. 

In general, we have to compute integrals over iV points, 
equivalent to NiN — 1) distances x, y, . . ., of the form 



1^ = 



inax(a;,y,...)<L 



. . .) 



(Jl) 



with a homogeneous function / such that the integral 
has a conformal weight (dimension in L) k: /(e) ^ L'^. 
For the integrals which appear in iV-loop diagrams, this 
weight is simply k = 2iVr (renormalization of g), or k 
2(A^ + \)t (renormalization of a). 

The integral over the distances is defined by the inte- 
gral over — 1 points, keeping one chosen point fixed. 
The residue is extracted from the dimensionless integral 



(J2) 



,f ■.= kL-'^1^ = L-'^L^l^ 
oL 

= L / f{x,y,...) max(x,y, . . .)" 

J max(a;,y,...)— L 



Note that the last factor is nothing but L~'^, but we have 
chosen this form as will become apparent shortly. The 
domain of integration can be decomposed into "sectors" , 
for instance 

{...<y<x = L},{...<x<y = L}, (J3) 

and we can map these different sectors onto each other 
by global conformal transformations. As we show below, 
we can for instance rewrite the integral (|J3I) sector by 
sector as 

J^eel/ /(x, ?/,...) max(a;,y, ...)-'' 

J x—L; y,... 

= L /(x, ?/,...) max(a;,y, ...)"'' . (J4) 

'J y—L; X,... 
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The constraint on the maximum of the distances is re- 
placed by the constraint on an arbitrarily chosen dis- 
tance. Due to our normalization f := tt- f d^z, this 
is equivalent to fixing the both endpoints of this largest 
distance. 

To formalize the above, consider the integral over a 
function / at order iV — 1 loops. Suppose f{zi, . . . , zn) 
is a homogeneous function of dimension — 2(7V — 1) -|- k. 
Define the function 



f{Zl,. . . ,ZN,) ■■= f{zi,...,ZN,) 



max{ I Zi 



;|} 



(J5) 

(This is the combination which appears in (|J2[) . Then 
the integral over zi, . . . , zj^-i (the relative coordinates 
between points), cut off by C(2i, . . . , zn) ■= Y[i<j 

z,\<L) 



In [a, L) 



f{zi,...,ZN)C{zi,...,ZN) (J6) 



Zl,...,Zjv_l 



has L-dimension 0. Consider a sector S (ordering of the 
distances). Be Xa \zi — Zj\, with 1 < a < m :— N{N — 
l)/2. Then S := {zi,...,}, s.t. xi < X2 < ■ ■ ■ < Xm- 
(Actually, we have chosen the labeling of the distances 
Xa to account for the ordering. This is not always the 
most practical thing to do.) Also define the character- 
istic function xsixi, . . . ,a;,„) of a sector S as being 1 if 
all distances satisfy the inequalities of the sector and 
otherwise. The L-derivative of the integral restricted to 
the sector S is 



J- 



L-^4{a,L) 



Xm—L 



(J7) 



The conformal mapping theore m^^i^^i^^ , whose proof we 
reproduce below for completeness, states that if the in- 
tegral (|J7I) is Riemann-integrable everywhere (or at least 
via finite-part prescription), then 



J'^ = I f(zi,---) xs{xi,.-.,Xm) 



(J8) 



In words: The above integral can be evaluated by fixing 
any of the distances to be L (or 1 equivalently). (|J4|) is 
then a simple corollary of (|J8|) . 



To prove the latter, we start from (|J7|) . First of all, 
since Xm = and introducing a (5-function to enforce it, 
becomes 



f{zi, . . . )6{Xm - L)XS{X1,- ■ ■,Xm) 



(J9) 



We now aim at integrating over the distances xi, . . . , Xm 
instead of the coordinates. For an arbitrary function g 
of the latter distances, this is 



d^zi . . . d^ZN^i gixi, . . . , Xm) 



d2;i...da;„ fi{xi, . . . , Xm)gixi, . . . , Xm) ■ (JIO) 
The measure is easily constructed as 

^(cci, . . . ,x™) = y"*^^^! ■ • -d^-^w-i S{xi - \zi - Z2I) X 

X . . .(5(a;„ - |zAr_i - ZAr|) , (Jll) 

where the (5-distributions enforce the x^'s to be the dis- 
tances between the Zj^s. 

We now want to map onto xi = L. To achieve this, we 
can always do the integration over a;; last. This gives for 

J'^ = J ^^i J dxi . . .dxi^idxi+i . . .dx„i fi{xi, . . . ,Xm) 

x6{Xjn~ L) f{xi,...,Xrn)Xs{xi,...,Xm) ■ (J12) 

We now make a change of variables. For all i but I, set 
Xi := XiXi/L . (J13) 



We also define xi :— L, and introduce this into (IJ12|) as 
I = J dxi 5{xi ~ L): 



J-^ = y dxi J dxi . . . dxjn ■ ■ • ,^m) S{xi - L) 

X /(ii, . . .,Xm)Xs{xi,. ■ .,Xm) 

X 5{imXi/L - L) — . (J14) 
xi 

Note that the factor of ^ consists of (77)^ 
from the terms dxi but di/; a factor of (77)^^ ^''^^ ^ 

from the measure; and a factor of (^) ^''^ 
Using that 



from /. 



dxi5{x„iXi/L- L)— , (J15) 
xi 



we obtain 

^ I dxi . . . dx„i fi{xi, . . . , x,n) S{xi - L) 

X /(ii, . . . , Xm)Q{Xm/xi < L / a) 
xxs{xi,. ■ . ,im) . (J16) 



Dropping the tildes, this is nothing but (|J12p with x„ 
replaced by xi which completes the proof. 
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